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ABSTRACT 


A  computer-aided  data  acquisition  system  was  developed 
and  a  microthermocouple  probe  constructed  to  obtain  thermal 
distributions  in  turbulent  buoyant  jets  exposed  to  a  cross- 
flowing  ambient  fluid.   The  system  performed  high  speed 
temperature  measurements  as  a  microthermocouple  probe  was 
automatically  traversed  through  a  sequence  of  preprogrammed 
positions  under  the  control  of  a  microcomputer.   Operability 
of  the  apparatus  was  demonstrated  by  measuring  temperature 
distributions  in  planes  perpendicular  to  the  streamwise  axis 
of  jets  from  which  contour  plots  of  temperature  were  gener- 
ated.  Using  temperature  distributions  along  with  velocity 
distributions  allow  buoyant  jet  characteristics  to  be 
computed,  including  the  entrainment  rate  of  ambient  fluid, 
jet  trajectory,  and  heat  transfer  to  the  ambient.   The 
experimental  technique  is  discussed  and  temperature  contour 
plots  for  a  jet  at  various  planes  are  presented. 
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NOMENCLATURE 


A.  .  Incremental  Cross-Sectional  Area  in  the 

-1  Temperature  Matrix 

B  Jet  Half-width 

b  Normalized  Jet  Half-width 

c  Specific  Heat 

P  ^ 

D  Diameter  of  the  Jet  at  the  Nozzle 

D.D  Binary  Mass  Diffusion  Coefficient 

F  Densiometric  Froude  Number 

g  Acceleration  of  Gravity 

Q  Heat  Transfer  Rate  from  the  Jet  to  the  Ambient 
Fluid 

R  Ambient- to-Nozzle  Flow  Ratio 

R=  Length  of  the  Probe  Arm 

a. 

r  Radial  Distance  from  the  Center  of  the  Jet 

r  Length  of  the  Probe 
P 

S  Schmidt  Number 

s  Streamwise  Coordinate  Along  the  Jet  Centerline 

T  Normalized  Jet  Temperature 

T  Ambient  Fluid  Temperature 

T. .  Jet  Temperature  Within  the  Temperature  Matrix 

T  Centerline  Jet  Temperature 

T  Nozzle  Temperature 

T  Jet  Temperature  As  Measured  by  the  Probe 

T(r)  Temperature  Within  the  Jet  at  a  Radial  Distance 
(r)  from  its  Center 


U  Centerline  Velocity  of  the  Jet  at  the  Nozzle 

U  Discharqe  Velocity  of  the  Jet 
o  r  j 

U(r)  Velocity  within  the  Jet  at  a  Radial  Distance  "r" 
from  its  Center 

U. .  Jet  Velocities  Corresponding  to  Locations  within 

1-'  the  Temperature  Matrix 

u  Normalized  Centerline  Velocity 

a  Entrainment  Coefficient;  Offset  Angle  of  the 
Probe  Arm 

3  Offset  Angle  of  the  Probe  Mounting  Bracket 

Y  Probe  Angle  of  Deflection  from  Horizontal 

0  Local  Angle  of  Inclination  from  Horizontal  of  the 
Jet  Streamwise  Axis 

A  Spreading  Ratio 

v  Kinematic  Viscosity 

p  Density  of  the  Jet  Fluid 

p  Density  of  the  Ambient  Fluid 
a 

p  Density  of  the  Ambient  Fluid  at  the  Nozzle  Exit 

Kao  ■* 

p  Density  of  the  Jet  at  Centerline 

m  2 

p  Centerline  Density  of  the  Jet  at  the  Nozzle 
'mo  J 

cj)  Angle  of  Inclination  of  the  Data  Plane  from  Horizontal 
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I.   INTRODUCTION 

Buoyant  jets  are  very  common  in  nature.   We  see  them  in 
the  form  of  exhaust  gases  emitted  from  smoke  stacks  of 
refineries,  mills  and  ships.   We  see  them  in  the  form  of 
heated  waste  water  expelled  into  the  sea  from  power  plants 
and  from  the  main  propulsion  condensers  in  steam  driven 
ships  and  submarines.   It  is  no  wonder  that  the  fluid 
mechanics  and  heat  transfer  characteristics  of  buoyant  jets 
have  been  of  interest  to  environmental,  civil  and  mechanical 
engineers  for  decades.   To  evaluate  their  ecological  impact, 
and  of  most  recent  interest,  to  harness  buoyant  jets  as  a 
means  of  detecting  military  targets  and  guiding  weapons,  it 
is  necessary  to  develop  models  which  accurately  predict 
their  trajectory  and  decay. 

Most  studies  to  date  have  dealt  with  buoyant  jets  rising 
through  a  quiescent  ambient  fluid;  however,  in  nature  most 
problems  involve  flowing  ambient  fluids.   Relatively  little 
experimental  work  has  been  done  with  buoyant  jets  in  cross- 
flow  and,  according  to  Hilder  [Ref.  1],  the  trajectories  of 
jets  and  the  entrainment  rates  of  ambient  fluid  predicted 
from  previous  work  do  not  agree  well  with  one  another.   Most 
mathematical  models  of  buoyant  jets  in  crossflow  assume 
Gaussian  profiles  for  velocity  and  temperature.   Nickodem 
[Ref.  2]  has  shown  through  experiments  that  in  fact,  the 
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Gaussian  profiles  of  velocity  are  altered  by  crossflow. 
This  leads  one  to  suspect  that  the  same  may  be  true  for  the 
temperature  profiles. 

The  objective  of  this  work  was  to  develop  a  system  to 
thermally  map  a  buoyant  jet  in  crossflow.   Then,  by  measur- 
ing both  velocity  and  temperature  distributions,  improved 
computations  of  entrainment,  trajectory  and  heat  transfer 
characteristics  of  jets  can  be  made  thereby  giving  rise  to 
more  accurate  models. 
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I I .   BUOYANT  JETS  DISCHARGED  TO  A  CROSSFLOW 

A.   PROPERTIES  QF  BUOYANT  JETS 

A  buoyant  jet  is  characterized  by  a  momentum  and  a 
density  differential  between  the  jet  and  its  surrounding 
ambient  resulting  from  a  variation  in  temperature  and/or 
fluid  concentrations.   Therefore,  fluid  motion  in  the  jet 
is  governed  by  both  inertial  and  buoyant  forces.   The  non- 
dimensional  ratio  of  these  forces,  known  as  the  densiometric 
Froude  number,  provides  an  important  quantitative  measure- 
ment of  jet  characteristics  and  is  shown  below. 

u2 

F   = 


gD(pa-po)/po 


where  U   is  the  jet's  discharge  velocity,  g  is  the  accelera- 
tion of  gravity,  D  is  the  discharge  diameter  of  the  jet, 

p   is  the  density  of  the  crossflowing  ambient  and  p_  is 
a  J  o 

the  density  of  the  jet  fluid  at  its  point  of  discharge. 

The  Gaussian  velocity  and  temperature  profiles  assumed 
by  most  models  of  buoyant  jets  are  very  similar.  Velocity 
behavior  is  given  by: 


U(r)   =   Um  exp(-r2/B  ) 


where  U   is  the  centerline  velocity,  r  is  the  independent 
m 

variable  and  a  radial  distance  from  the  centerline  of  the 
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jet,  and  B  is  defined  as  nominal  jet  halfwidth.   As  r 

approaches  B,  velocity  decays  to  (l/e)U   [Ref.  3].   Simi- 

m 

larly,  temperature  behavior  is  given  by: 


T(r)   =   Tm  exp(-r2/A2B2: 


where  T   is  the  centerline  temperature,  r  and  B  are  defined 
m 

the  same  as  above  and  A,  a  spreading  ratio,  is  the  inverse 
of  the  turbulent  Schmidt  number  (S) .   S  is  defined  as  the 
ratio  of  the  molecular  momentum  and  mass  dif f usivities  and 
is  equal  to  v/D   where  v  is  the  kinematic  viscosity  and 
D__,  is  the  binary  mass  diffusion  coefficient  associated  with 
substances  A  and  B  [Ref.  4 J .   Although  A  varies  inversely 
with  the  Froude  number,  the  change  is  very  slight,  and  in 
the  case  where  substances  A  and  B  are  both  water,  A  is 
slightly  greater  than  1.   Hirst  [Ref.  3]  found  A  to  vary 
between  1.16  at  F  =  0  to  1.11  at  F  =  infinity.   The  net 
effect  then,  is  a  more  gradual  temperature  decay  than  was 
found  with  velocity. 

Most  buoyant  jet  models  consider  the  entrainment  of  the 
ambient  fluid  into  the  jet  and  are  based  on  relevant  conser- 
vation equations  of  mass,  momentum  and  energy.   In  conserva- 
tion of  mass,  the  downstream  change  in  total  mass  of  the  jet 
is  equated  to  the  mass  of  the  entrained  fluid.   The  conser- 
vation of  momentum  must  consider  both  vertical  and  horizontal 
contributions.   Changes  in  vertical  momentum  are  equated 
to  the  buoyant  forces  while  changes  in  the  horizontal 
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momentum  of  the  jet  are  equated  to  the  horizontal  momentum 
of  the  entrained  fluid.   The  conservation  of  energy  involves 
energy  changes  resulting  from  variations  in  the  ambient 
temperature  as  caused  by  the  jet.   Hilder  [Ref.  1]  developed 
the  following  governing  equations  in  non-dimensional  differ- 
ential form. 


CONTINUITY  ^-(u^b2)   =   2ab  [  |u  -R  cos  8  |  +  a^sin  9] 


HORIZONTAL  MOMENTUM     ^-(u2b2cos8)   =   4Rab [ I u  -Rcosel  +a,Rsinf 

ds   m  '  m       '3 


VERTICAL  MOMENTUM       ^-(u2b2  sin8)  =   ( — )  ■  2X   i? 

ds   m  p   -p      _2 

Kao   mo    F 

-,                   ,2  2  At      0 

ENERGY                    2_(u  T'-Kp- )  =   -r^(u  b^} 

ds   m   ,  ,  2   .  ds   m 
{  A  +1; 


B.   FLOW  REGIMES 

The  jet  passes  through  several  regimes  as  it  travels 
from  the  nozzle  through  the  ambient.   The  three  regions  most 
frequently  referred  to  are  shown  in  Figure  1.   They  are  the 
zone  of  flow  establishment,  the  zone  of  established  flow 
and  the  far-field  zone  [Ref.  3].   In  the  zone  of  flow  estab- 
lishment, the  velocity  and  turbulence  profiles  transform  from 
the  conditions  within  the  nozzle  to  a  free  turbulent  flow 
condition.   It  is  in  this  region  that  the  jet  begins  to  mix 
with  the  ambient  fluid;  however,  the  flow  is  still  more 
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strongly  influenced  by  the  nozzle  discharge  conditions  than 
by  the  ambient.   When  the  turbulent  mixing  has  reached  the 
centerline  of  the  jet,  the  zone  of  established  flow  is  said 
to  begin.   In  this  region,  the  profiles  have  assumed  their 
free  turbulent  shapes.   Now  the  flow  is  governed  by  the 
jets'  momentum  and  buoyancy  as  well  as  by  the  condition  of 
the  crossflow.   The  far  field  zone  is  defined  as  that  region 
in  which  jet  momentum  is  depleted  and  the  jet  fluid  is  con- 
vected  and  diffused  by  the  ambient  currents  and  turbulence. 

C.   EFFECTS  OF  CROSSFLOW 

At  the  immediate  exit  of  a  cylindrical  nozzle,  a  verti- 
cally discharged  buoyant  jet  has  a  nearly  uniform  velocity 
distribution  and  has  the  same  cross-sectional  shape  as  the 
nozzle  itself.   The  velocity  gradient  between  the  jet  and 
the  crossflowing  ambient  creates  longitudinal  shear  stresses 
at  the  jet's  sides,  a  positive  pressure  region  immediately 
upstream  and  a  negative  pressure  region  immediately  downstream 
of  the  jet.   This  results  in  the  deflection  of  the  jet's 
trajectory  in  the  downstream  direction  (Figure  2),  the 
creation  of  counterrotating  vortices  at  the  jet's  outer 
edges  and  the  deformation  of  the  original  circular  cross- 
sectional  shape  into  the  form  of  a  kidney.   As  the  stream- 
wise  axis  of  the  jet  approaches  the  direction  of  the  cross- 
flow,  these  effects  become  progressively  less  pronounced. 
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III.       EXPERIMENTAL    APPARATUS 

A.  SYSTEM  OVERVIEW 

A  surplus  milling  machine  was  configured  with  synchron- 
ous drive  motors  interfaced  with  a  microcomputer  that  auto- 
matically positioned  its  bed.   It  was  used  as  a  three-dimensional 
positioning  platform  in  the  same  manner  as  in  the  laser 
Doppler  velocimetry  work  undertaken  by  Nickodem  [Ref.  2]. 
The  milling  machine  was  placed  adjacent  to  a  rectangular 
plexiglass  flume  through  which  the  crossflowing  ambient 
fluid  flowed.   A  vertical  nozzle  was  installed  in  the  base 
of  the  flume  to  provide  the  jet.   A  temperature  probe  was 
suspended  through  an  opening  in  the  top  of  the  flume  above 
the  nozzle  by  an  arm  attached  to  a  base  mounted  on  the 
milling  machine  bed  as  shown  in  Figure  3.   As  the  probe  was 
automatically  traversed  through  a  series  of  preprogrammed 
positions  across  the  jet,  temperature  data  was  automatically 
sensed  and  stored  at  high  speeds  by  the  computer. 

B.  CROSSFLOW  SYSTEM 

As  illustrated  in  Figure  4,  the  crossflow  circulation 
pump  took  water  from  the  cylindrical  24  8.8  1  (65.7  gal) 
reservoir  shown  in  Figure  5  and  discharged  through  5.076  cm 
(2  in)  diameter  tubing  into  a  cylindrical  flow  settling 
chamber  30.46  cm  (12  in)  in  diameter  and  60.91  cm  (24  in) 
tall  located  within  the  60.91  cm  x 60. 91  cm  *  88.83  cm 
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(24  in  x 24  in  *  35  in)  inlet  chamber  shown  in  Figure  6.   The 
settling  chamber  was  sealed  at  its  bottom  so  that  the  water 
spilled  from  its  top  into  the  inlet  chamber  through  honey- 
combed flow  straighteners  to  reduce  turbulence  and  evenly 
disperse  the  flow.   To  further  reduce  turbulence,  the  flow 
was  broken  by  another  stack  of  honeycombed  flow  straighteners 
and  a  layer  of  fiberglass  filter  material  located  immediately 
above  the  normal  operating  water  level.   The  flow  next 
entered  a  24.4  cmx32.39  cm  xl82.9  cm  (9.625  in  xl2.75  in  x  72  in! 
flume  shown  in  Figure  7  through  a  vertical  section  of  the 
same  honeycombed  material  mentioned  above.   To  avoid  inadver- 
tent spillage  over  the  sides  of  the  flume  during  system 
start-up,  a  5.076  cm  (2  in)  diameter  overflow  pipe  was 
located  in  the  inlet  chamber.   During  normal  operation, 
a  gate  valve  in  this  piping  was  closed.   The  flow  left  the 
flume  through  a  7.614  cm  (3  in)  diameter  pipe  at  its  end 
and  re-entered  the  crossflow  circulation  pump  reservoir.   A 
gate  valve  located  in  this  piping  and  shown  in  Figure  8  was 
used  to  regulate  the  water  level  and  flow  velocity  in  the 
flume.   The  optimum  adjustment  of  this  valve  was  determined 
by  trial  and  error  to  be  closed  two  turns  from  its  fully 
open  position.   Either  a  globe  or  ball  valve  would  have  been 
more  appropriate  for  this  purpose;  however,  neither  was 
readily  available,  so  the  gate  valve  was  used.   The  bracket 
shown  at  the  base  of  the  flume  in  Figure  8  maintained  align- 
ment between  the  flume  and  the  milling  machine.   The  water 
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Figure  7.   Flume  Arrangement 


Figure  8.   Flume  Discharge  Piping 
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in  the  flume  was  gradually  heated  by  repetitive  circulation 
through  the  crossflow  pump  and  by  the  addition  of  the  heated 
water  from  the  jet.   To  maintain  a  constant  temperature 
crossflow,  cooling  water  from  a  refrigerated  bath  shown  in 
Figure  9  was  circulated  through  a  coil  of  1.269  cm  (.5  in) 
diameter  copper  tubing  located  in  the  flow  settling  chamber. 
Also,  fresh  water  was  added  at  the  flow  settling  chamber 
as  an  equal  amount  was  drained  from  the  crossflow  pump 
reservoir  through  a  1.269  cm  (.5  in)  diameter  pipe.   Cross- 
flow  temperature  was  monitored  by  a  Type-T  thermocouple 
located  in  the  inlet  chamber. 

C.   JET  SYSTEM 

In  reference  to  Figures  4  and  5,  the  jet  flow  pump 
circulated  water  from  a  rectangular  26.27  1  (6.94  gal) 
reservoir  and  discharged  through  1.26  cm  (.5  in)  diameter 
tubing  to  a  33.0  cm  x  50 . 8  cmx54.6  cm  (13  in  x20  in  x  21 . 5  in) 
head  tank  (Figure  10) .   The  amount  of  flow  to  the  head  tank 
was  regulated  by  a  globe  valve.   Due  to  a  low  flow  rate 
to  the  head  tank,  water  was  also  recirculated  back  to  the 
reservoir  in  order  to  maintain  sufficient  flow  through  the 
jet  pump  to  prevent  overheating  it.   A  constant  water  level 
was  maintained  in  the  head  tank  by  a  stand  pipe  which  allowed 
overflow  back  to  the  reservoir.   Sufficient  flow  into  the 
tank  was  maintained  to  make  sure  that  it  slightly  overflowed 
continuously.   Water  drained  from  the  bottom  of  the  head 
tank  through  1.26  cm  (.5  in)  diameter  tubing  and  passed 
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Figure  9.   Refrigerated  Bath 
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Figure  10 .   Head  Tank 
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through  a  rotometer,  a  .95  cm  (.375  in)  tubing  reducer,  a 
dye  injection  system  (Figure  11) ,  a  water  heater  (Figure  12) 
consisting  of  approximately  6.09  m  (20  ft)  of  .95  cm  (.375 
in)  diameter  copper  tubing  coiled  in  a  heated  bath  and 
finally  a  7.144  cm  (.28125  in)  nozzle  which  discharged  into 
the  bottom  of  the  flume.   Flow  was  controlled  by  pinching 
the  tubing  between  the  heater  and  the  nozzle  with  surgical 
clamps.   Drainage  from  the  crossflow  reservoir  discussed  in 
Section  III.B  was  used  to  replenish  the  jet  reservoir.   The 
dye  injection  system,  used  in  photographing  the  jet,  was 
located  approximately  8.23  m  (27  ft)  upstream  of  the  nozzle 
to  minimize  any  disturbance  to  the  jet  that  it  might  have 
caused.   The  majority  of  this  distance  was  taken  up  by  the 
heating  coil  mentioned  above.   The  vertical  distance  between 
the  top  of  the  stand  pipe  in  the  head  tank  and  the  tip  of 
the  nozzle  in  the  flume  was  2.2  m  (86.5  in)  which  equated  to 
21.56  KPa  (3.127  psig) .   Jet  flow  temperature  was  monitored 
by  a  Type-T  thermocouple  located  within  the  jet  flow  tubing 
approximately  1.167  m  (46  in)  from  the  nozzle. 

D.   TEMPERATURE  PROBE 

Measuring  temperatures  in  a  buoyant  jet  with  a  thermo- 
couple is  intrusive.   To  reduce  the  probability  of  distorting 
results,  steps  were  taken  to  minimize  the  cross-sectional 
area  of  the  temperature  measuring  device  as  seen  by  the  flow 
of  the  jet.   A  .0254  mm  (.001  in)  diameter  Type-E  micro- 
thermocouple  was  selected.   The  suspension  device  for  the 
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microthermocouple  had  to  be  rigid  and  have  a  small  cross- 
sectional  area,  for  reasons  discussed  above,  as  well  as  be 
an  electrical  insulator  to  prevent  interference  with  the 
thermocouple  performance.   A  glass  annulus  approximately 
1.45  mm  (.057  in)  in  diameter  and  11.27  cm  (4.4375  in)  in 
length  was  chosen.   One  lead  of  the  thermocouple  was 
threaded  through  the  annulus  and  the  other  was  glued  with  a 
fast  drying  modelers'  glue  along  the  outer  surface,  allowing 
the  microthermocouple  junction  to  protrude  slightly  from  the 
tip  of  the  annulus.   The  leads  at  the  opposite  end  of  the 
annulus  were  welded  to  .0762  mm  (.003  in)  diameter  wire 
which  subsequently  was  connected  to  28  AWG  extension  wire 
to  the  computer.   The  annulus  was  mounted  as  shown  in  Figure 
13.   Henceforth,  this  device  will  be  referred  to  as  the 
probe . 

E.   PROBE  ACTUATOR-  ASSEMBLY 

The  cross  sectional  area  of  the  probe  as  seen  by  the  jet 
was  further  reduced  by  orienting  the  probe  tangentially  to 
the  trajectory  of  the  jet  as  shown  in  Figure  14.   This 
photograph  indicated  that  the  probe  created  no  noticeable 
interference  with  the  jet  hydrodynamics.   Probe  orientation 
was  accomplished  by  the  linkage  assembly  shown  in  Figure  15. 
The  fixed  end  of  the  probe  was  hinged  to  a  streamlined  tube 
23.495  cm  (9.25  in)  long  with  a  maximum  width  and  depth,  as 
seen  by  the  jet,  of  3.175  mm  (.125  in)  and  6.35  mm  (.25  in) 
respectively.   It  was  rigidly  connected  to  the  mounting 
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Figure  13.   Probe  Profile 


Figure  14.   Probe  in  Jet 
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A.  Drive  Motor  Pulley 

B.  Cable  and  Turnbuckle 

C.  Stroke  Pulley 

D.  Linkage  Rod 

E.  Mounting  Base 

F.  Probe  Arm 

G.  Mounting  Bracket 
H.  1/2-13  Nut 
I.  1/2-13  Hollowed  Stud 
J.  Streamlined  Tube 
K.  Hinge  Assembly 
L.  Probe   (Glass  Annulus) 


Figure    15.       Probe   Assembly   Linkage 
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bracket  as  shown  in  Figure  15  which  was  connected  to  the 
probe  arm  shown  in  Figure  3  by  a  single  stud  which  allowed 
pivoting  of  the  probe  from  side- to-side .   The  stud  was  also 
hollowed  so  that  a  linkage  rod  could  extend  from  the  hinge 
assembly  through  the  tube  and  stud  to  a  stroking  pulley 
which  was  rotated  by  a  small  motor.   The  hinge  assembly  and 
the  stroking  pulley  were  spring  loaded  to  reduce  hysteresis. 
As  shown  in  Figures  16  and  17,  the  1.5  VDC  motor,  geared  to 
one  rpm,  was  directly  coupled  to  a  potentiometer  as  well  as 
the  drive  pulley.   The  potentiometer  was  configured  in  a 
voltage  divider  such  that  the  amount  of  motor  rotation,  and 
ultimately  the  degree  of  probe  deflection,  was  proportional 
to  the  potential  difference  sensed  across  the  potentiometer. 
Limit  switches  were  installed  at  the  stroke  pulley  as  shown 
in  Figures  16  and  18  to  prevent  damage  to  the  linkage 
assembly  due  to  over-rotation. 

F.   MICROCOMPUTER  INTERFACE 

The  data  collection  process  consisted  of  adjusting  the 
probe  angle  of  deflection,  traversing  the  three-dimensional 
positioning  platform  and  measuring  temperature  profiles. 
All  of  the  mechanisms  which  controlled  these  events  were 
interfaced  to  an  HP-9  826  computer  shown  in  Figure  19  through 
an  HP-6942A  multiprogrammer  which  performed  high  speed  analog- 
to-digital  conversions  and  ultimately  provided  control 
signals  to  govern  relays  within  the  system.   Refer  to  MAIN_T 
in  Appendix  B  for  the  microcomputer  software  which  directed 
this  process . 
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Figure  17.   Probe  Actuator  Motor- 
Potentiometer  Arrangement 


Figure  18.   Probe  Actuator 


35 


Figure  19.   HP-9826  Microcomputer 
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1 .  Probe  Angle  Adjustment 

As  discussed  in  Section  III.E,  a  potentiometer  con- 
figured as  a  voltage  divider  provided  probe  angle  feedback 
to  the  computer  as  shown  in  Figure  20.   The  direction  of 
motor  rotation  was  controlled  by  the  computer  through  Type- 
MRR1CDL  replays  connected  as  shown  in  Figure  21.   When  the 
probe  was  at  a  desired  position,  a  5.0  VDC  signal  was  applied 
to  pins  6L  and  6R  which  allowed  both  relays  to  assume  the 
normally  closed  (NC)  positions  which  opened  the  power  circuit 
to  the  motor.   When  it  was  desired  to  rotate  the  motor  clock- 
wise, pin  6R  was  grounded  which  resulted  in  5.0  VDC  applied 
across  the  coil  in  the  right-hand  relay.   This  caused  the 
relay  to  assume  its  normally  open  (NO)  position  resulting  in 
a  1.5  VDC  signal  at  terminal  B  of  the  motor  causing  it  to 
rotate  in  the  clockwise  direction.   The  left-hand  relay  was 
activated  in  a  similar  manner  for  counterclockwise  rotation. 
Refer  to  PROBE_SUBS  in  Appendix  B  for  probe  positioning 
software . 

2.  3-D  Positioning  Platform  Movement 

Positioning  platform  movement  was  controlled  in  a 
manner  similar  to  the  probe  and  was  discussed  in  detail  by 
Nickodem  [Ref .  2] .   Refer  to  MTR_SUBS  in  Appendix  B  for 
associated  HP-9826  software. 

3 .  Temperature  Data  Collection 

Three  thermocouples  were  monitored  in  the  data  col- 
lection process.   A  Type-T  thermocouple  located  in  the  inlet 
chamber  measured  the  ambient  fluid  temperature  in  the  flume, 
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a  Type-T  thermcouple  located  in  the  tubing  between  the  heater 
and  the  nozzle  measured  nozzle  temperature  and  a  Type-E 
microthermocouple  in  the  probe  measured  the  temperature  in 
the  jet.   The  EMF's  generated  by  these  thermocouples  were 
amplified  by  "Omega  Omni-Amp  IIB"  millivolt  amplifiers  shown 
in  Figure  22  prior  to  entering  the  multiprogrammer  for 
analog-to-digital  conversion  and  eventual  transformation  to 
temperature  readings.   Fourth-order  least  squares  coeffi- 
cients for  this  conversion  were  taken  from  Beckwith  [Ref .  5] . 
Operation  of  the  crossflow  circulation  pump  created  suffi- 
cient electrical  interference  to  distort  the  thermocouple 
signals.   This  problem  was  corrected  by  applying  a  thin 
coating  of  silicon  sealant  to  the  Type-T  thermocouple  junc- 
tions and  by  connecting  the  crossflow  circulation  pump 
casing,  the  nozzle  and  the  jet  tubing  in  the  vicinity  of  the 
nozzle  thermocouple  to  a  common  ground.   Because  the  jet 
tubing  was  plastic,  it  was  necessary  to  manufacture  a  brass 
"T"  connector  as  shown  in  Figure  23  which  was  grounded  and 
located  in  close  proximity  to  the  thermocouple  junction. 
Refer  to  T  SUBS  in  Appendix  B  for  associated  software. 
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Figure  22.   Thermocouple  Amplifiers  and  HP-6942A 
Multiprogrammer 
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Figure  23.   Nozzle  Thermocouple  Grounding  Arrangement 
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IV.   EXPERIMENTAL  PROCEDURES 

A.  CALIBRATION 

Before  the  data  collection  process  could  begin,  the 
rotometer,  the  thermocouples,  the  probe  and  the  positioning 
platform  had  to  be  calibrated.   The  detailed  steps  taken 
are  discussed  below. 

1 .  Rotometer 

With  a  constant  level  maintained  in  the  head  tank 
and  the  jet  tubing  disconnected  from  the  nozzle  and  elevated 
to  the  same  height  as  the  top  of  the  nozzle,  five  100  ml 
samples  were  drawn  through  the  rotometer  and  timed  to  the 
nearest  0.01  second  at  each  rotometer  reading  from  10% 
to  75%  in  5%  increments.   Flowrates  and  standard  deviations 
in  ml/s  are  shown  in  Table  1. 

2 .  Thermocouples 

Since  nozzle  and  probe  temperatures  were  to  be 
normalized  by  the  ambient  temperature,  the  nozzle  and  probe 
thermocouples  were  calibrated  relative  to  the  ambient  thermo- 
couple by  using  the  microcomputer  program  T_CAL  in  Appendix 

B.  The  procedure  followed  is  outlined  in  the  initial  com- 
ments of  the  program.   Coefficients  for  first  order  curve 
fits  were  solved  by  the  least  squares  method  with  the 
mainframe  programs  TCAL  and  TFIT  found  in  Appendix  C. 
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3 .  Probe 

As  the  probe  assembly  was  being  developed,  it  was 
convenient  to  test  its  suitability  with  the  probe  calibration 
panel  shown  in  Figure  24.   Resistance  changes  across  the 
potentiometer  were  recorded  for  varying  degrees  of  deflec- 
tion.  Analysis  of  this  information  led  to  improved  designs 
from  the  standpoint  of  reduced  hysteresis  and  repeatability. 
The  microcomputer  program  PROBE_CAL  in  Appendix  B  was 
developed  to  enable  calibration  of  the  final  design  after 
it  was  installed  in  the  system  as  shown  in  Figures  25  and 
26.   The  calibration  procedure  is  outlined  in  the  preliminary 
comments  of  the  program. 

4 .  3-D  Positioning  Platform 

The  positioning  platform  was  calibrated  in  a  manner 
that  placed  the  tip  of  the  probe  at  desired  locations  within 
the  flume  relative  to  the  tip  of  the  nozzle.   Referring  to 
the  coordinate  system  illustrated  in  Figure  2,  the  center 
of  the  nozzle  was  defined  as  (0,0,0)  in  xyz-coordinates . 
The  following  relationships  apply  to  the  probe  geometry 
shown  in  Figure  27 : 


X(real)       =      X      -    R     cos  a    tan(iT/4    -    a/2) 
o  a 


Y(real)       =      Y      +   r    (1   -    cos    y)    +    R      cos    a 
op  a 


Z(real)       =      Z      -    r      sin   y 
o  p  ' 
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Figure  26.   Probe  Calibration  Panel 
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SIDE  VIEW 


Flume 


Nozzle 


Probe 


TOP  VIEW 


(h- *-- 
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Probe  Mounting 
Bracket 


(1/8  scale) 


Figure  27.   Probe  Positioning  Geometry 
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(X  ,Y  ,  Z  )  was  the  position  of  the  milling  machine  bed 
o   o   o  ^  r 

when  the  tip  of  the  probe  was  at  position  (0,0,0)  with 

Y  =  0  and  a  =  6  =  90  degrees.   For  calibration,  the  probe 

and  probe  arm  were  configured  with  these  settings  as  shown 

in  Figure  28.   With  measured  values  of  r   and  R   entered 

*  pa 

into  the  microcomputer  program  MAIN_T,  the  calibration  was 
accomplished  by  the  program  MOTOR_CAL.   The  step-by- step 
procedure  followed  was  outlined  in  the  subprogram  "SUB 
Calibrate."   As  illustrated  in  Figures  3,  27  and  28,  the 
length  of  R   could  be  modified  to  compensate  for  adjustments 
of  a  and  3  to  positions  other  than  90  degrees.   Decreasing 
a    increased  the  distance  along  the  Y-axis  in  which  the 
probe  could  be  positioned.   Increments  of  a  and  (3  were 
scribed  on  the  top  of  the  base  and  at  the  tip  of  the  probe 
arm  in  Figure  2  8  to  accommodate  this  change,  if  desired. 
The  program  MAIN_T  queried  the  user  for  the  value  of  a  and 
assumed  a   =  £.   The  calibration  software  also  established 
position  limits  to  prevent  driving  the  probe  into  the  sides 
of  the  flume. 

B.   PRELIMINARIES 

Crossflow  velocity  was  determined  by  injecting  blue  food 
coloring  into  the  flow  and  timing  its  travel  through  a  1.0  m 
interval.   The  average  of  several  trials  indicated  the 
velocity  was  .130  m/s  (.427  ft/sec)  with  the  flume  outlet 
valve  closed  two  turns  from  its  fully  open  position. 
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Figure  28.   Probe  Arm  Positioning  for  Motor  Calibration 
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Figure  29.   Typical  Buoyant  Jet  as  Obsered  with  Dye 
Injected 
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Photographs  of  jet  profiles  as  shown  in  Figure  29  were 
taken  to  determine  jet  trajectory  and  halfwidths  along  the 
streamwise  axis.   This  was  done  by  injecting  blue  food 
coloring  into  the  jet  flow  as  discussed  in  Section  III.C. 
The  specific  gravity  of  the  food  coloring  was  found  to  be 
considerably  less  than  that  of  water.   To  eliminate  the 
added  buoyant  effect  this  would  have  had  on  the  jet,  a  small 
quantity  of  alcohol  was  mixed  with  the  food  coloring  as 
suggested  by  Merzkirch  [Ref .  6] .   The  amount  of  alcohol  added 
was  determined  by  trial  and  error.   As  small  quantities  were 
added  and  mixed,  samples  were  gently  placed  on  the  surface 
of  a  beaker  of  water.   Pure  food  coloring  laid  on  the  sur- 
face and  very  slowly  mixed  with  the  water.   As  alcohol  was 
added,  this  buoyant  effect  grew  progressively  less  and  the 
mixture  would  settle  into  the  water.   The  mixture  was  con- 
sidered satisfactory  when  it  no  longer  laid  on  the  surface, 
but  settled  to  some  equilibrium  position  in  the  beaker. 

Slide  photographs  of  the  jets  were  projected  onto  large 
sheets  of  3.175  mm  (.125  in)  grid  graph  paper  and  digitized 
along  approximate  streamwise  axes  and  half-width  trajectories 
A  scaling  factor  was  determined  by  equating  the  projected 
width  of  the  nozzle  to  its  known  outer  diameter  of  7.9375  mm 
(.3125  in).   The  above  data  was  fit  to  the  Michaelis-Menter 
Equation  [Ref.  7]  shown  below  by  the  least  squares  method 
with  the  mainframe  program  JETCURV  in  Appendix  C: 


b  +  y 
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Correlation  coefficients  close  to  1  were  consistently 
obtained.   To  determine  positions  within  the  jet  at  which 
to  make  temperature  measurements,  five  evenly  spaced  posi- 
tions per  jet  flow  rate  were  selected  along  the  streamwise 
axis  in  the  zone  of  established  flow.   Data  planes  slightly 
larger  than  the  jet  width  were  centered  at  these  points  and 
oriented  perpendicular  to  the  streamwise  axis.   One  hundred 
data  points  were  selected  in  a  symmetric  square  matrix  with 
points  most  densely  populated  near  the  center.   The  planes 
were  identified  alphabetically  and  in  consecutive  order  from 
"A"  to  "F" ,  where  "A"  represented  the  plane  nearest  the 
nozzle.   The  positions  were  entered  into  the  microcomputer 
and  stored  by  plane  on  a  floppy  disk  by  the  program  LOAD_XYZ 
in  Appendix  B.   Accompanying  each  data  point  was  a  probe 
deflection  angle  used  to  orient  the  probe  parallel  to  the 
path  of  the  jet  to  minimize  interference.   This  angle  was 
determined  by  evaluating  the  first  derivatives  of  the 
equations  developed  for  the  streamwise  axis  and  the  half- 
widths  and  performing  an  interpolation  based  upon  the  data 
points'  position  relative  to  the  two  curves. 

C.   DATA  ACQUISITION 

The  flow  systems  were  placed  into  operation  and  the 
ambient  and  nozzle  temperatures  were  monitored  with  the 
program  T_SUBS  to  evaluate  system  stability  and  readiness 
for  data  acquisition.   The  system  usually  took  approximately 
two  hours  to  come  into  equilibrium.   This  could  be  monitored 
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by  watching  the  jet  nozzle  temperature.   When  conditions 
were  stable,  data  consisting  of  two  hundred  probe,  ten 
ambient  and  ten  nozzle  temperature  samples  per  position  was 
collected,  one  plane  at  a  time,  by  the  program  MAIN_T. 
The  following  information  was  stored  on  a  floppy  disk  for 
each  data  point:   x,  y,  and  z  coordinates;  mean  probe, 
nozzle  and  ambient  temperatures  and  the  standard  deviation 
of  the  probe  measurements.   The  data  was  transferred  to 
the  mainframe  computer  by  using  a  modem,  the  microcomputer 
program  SEND_DATA  and  the  mainframe  program  GRAB. 

D.   DATA  REDUCTION 

The  raw  data  was  organized  into  a  more  usable  format 
by  the  mainframe  program  TDATA  which  also  converted  the  XYZ 
coordinates  into  the  XSW  system  shown  in  Figure  2.   The 
resulting  data  was  selectively  sent  to  the  program  C0NT0UR4 
which  applied  calibration  coefficients  to  the  temperature 
data  and  normalized  it  in  the  following  manner: 

T   -  T 

T   =  _P a 

T   -  T 
n     a 

where  T  was  the  jet  temperature  as  measured  by  the  probe, 

T  was  the  ambient  fluid  temperature  and  T  was  the 
a  ^  n 

temperature  of  the  jet  within  the  nozzle. 

Contour  plots  of  this  information,  generated  by  the 
CONTOUR  option  of  the  graphics  package  DISSPLA  [Ref.  8], 
are  presented  in  Figures  30  through  35. 
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Figure  30.   Plane  A  Temperature  Contour  Plot 
(large  scale) 
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Figure  31.   Plane  A  Temperature  Contour  Plot 
(small  scale) 
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Figure  32.   Plane  B  Temperature  Contour  Plot 
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TEMPERATURE  CONTOURS  IN  A  BUOYANT  JET 
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Figure 


33.   Plane  C  Temperature  Contour  Plot 
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Wm  k  CROSSFLOWSNQ  AMBIENT 
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Fiqure  34.   Plane  D  Temperature  Contour  Plot 
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Figure  35.   Plane  E  Temperature  Contour  Plot 
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E.   RESULTS 

The  contour  plots  found  in  Figures  30  through  35  support 
the  conjectures  in  Section  II. C  concerning  the  effects  of 
the  crossflow.   Figure  30,  a  plot  of  plane  A  centered  on 
the  streamwise  axis  46  degrees  from  horizontal  and  7.327  mm 
(.288  in)  downstream  of  the  nozzle  shows  significant  distor- 
tion.  Figures  31  through  35  show  planes  A  through  E  sequen- 
tially plotted  on  the  same  scale  in  order  to  observe  overall 
jet  behavior.   Plane  E  was  located  87.313  mm  (3.438  in) 
downstream  of  the  nozzle  and  86  degrees  from  horizontal. 
It  can  be  seen  that  as  the  jet  traveled  further  downstream, 
the  distorting  effect  grew  progressively  less,  as  expected. 

The  rate  of  heat  transfer  from  the  jet  to  the  ambient 
was  calculated  for  each  plane  utilizing  the  temperature 
distribution  matrix  generated  in  the  program  CONTOUR4  in 
Appendix  C  and  the  following  relationship: 


m    n 

q   =     y     y    p  c   a. .  u. .  t. . 
±£i  j=i      p    ^    ^    ^ 


where  p  was  the  relative  density  of  the  jet,  A. .  was  the  area 

of  each  matrix  segment,  U. .  was  the  velocity  in  each  segment, 

c   was  the  specific  heat  of  the  jet  and  T . .  was  the  tempera- 
P  i-J 

ture  in  each  segment.   Velocity  was  measured  along  the 
streamwise  axis  by  a  laser  Doppler  velocimeter  (LDV) .   It 
appeared  to  be  nearly  constant  in  the  region  of  the  jet 
observed.   The  mean  and  standard  deviation  was  44.875  mm/s 
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and  1.8  mm/s  respectively.   By  assuming  the  Gaussian  profile 
shown  in  Section  II. B,  velocity  was  determined  at  radii 
corresponding  to  each  segment  in  the  matrix  mentioned 
above.   The  rates  of  heat  transfer  are  shown  in  Table  2. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  objective  of  this  thesis  was  to  develop  a  computer- 
aided  data  acquisition  system  and  construct  a  microthermo- 
couple  probe  to  be  used  by  follow-on  students  to  study 
temperature  distributions  in  turbulent  buoyant  jets.   Sample 
data  was  taken  to  verify  system  operability.   Based  on 
results,  the  system  performed  in  a  satisfactory  manner  and 
will  be  an  invaluable  tool  for  subsequent  studies  of  buoyant 
jets  in  a  crossflowing  ambient. 

Data  point  positions  were  hand  calculated  and  loaded 
into  the  microcomputer  with  the  program  LOAD_XYZ.   This  was 
an  extremely  time  consuming  task  and  distracted  the  user 
from  defining  more  than  100  data  points  per  plane.   The 
system  can  be  greatly  improved  with  the  addition  of  micro- 
computer software  that  would  automatically  determine  and 
load  data  point  positions.   Data  point  population  could  then 
be  increased  with  ease  which  should  result  in  smoother  con- 
tour plots  and  more  accurate  heat  transfer  calculations. 

It  was  necessary  to  continuously  add  fresh  water  to  the 
crossflow  system  as  an  equal  amount  was  drained  from  it  in 
order  to  maintain  the  crossflowing  ambient  at  a  constant 
temperature.   This  was  because  the  cooling  coil  located 
within  the  flow  settling  chamber  was  inadequate  to  compen- 
sate for  the  heat  added  to  the  system  by  the  jet  and  the 


61 


crossflow  circulation  pump  by  itself.   Although  it  was 
possible  to  maintain  the  ambient  constant  within  1.4  C  by 
this  method,  in  the  interest  of  conserving  water,  it  is 
recommended  that  either  a  larger  capacity  chilled  water 
bath  or  a  cooling  system  that  circulates  a  refrigerant 
rather  than  water  be  appropriated  for  this  purpose. 
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APPENDIX  A 
UNCERTAINTY  ANALYSIS 

Experimental  uncertainty  was  analyzed  in  accordance  with 
the  guidelines  set  forth  by  Holman  [Ref .  9] .   Uncertainties 
in  the  primary  measurements,  based  upon  manufacturer  speci- 
fications and/or  the  number  of  significant  digits  which 
could  be  read,  follow: 

1.  Time  (rotometer  calibration):  .01  s 

2.  Volume  (rotometer  calibration):  1.0  ml 

3.  Rotometer  reading:  1.0  percent 

4.  3-D  positioning  platform  resolution:  .1524  mm 

5.  Probe  deflection  angle  resolution:  .5  degrees 

6.  Thermocouple  resolution:  .0099  °C 

7.  Thermocouple  time  constant:  .004  s 

8.  Time  (crossflow  velocity):  .01  s 

9.  Length  (crossflow  velocity) :  .05  m 

Based  on  the  above,  the  uncertainty  of  the  flow  rate  of  the 
jet  was  estimated  to  be  1.0  ml/s .   Maximum  uncertainty  in 
the  position  of  the  probe's  tip  due  to  the  uncertainty  in 
deflection  angle  was  determined  to  be  .983  mm  (.0387  in). 
Combined  with  the  resolution  of  the  positioning  platform, 
the  tip  of  the  probe  was  positioned  with  an  uncertainty  of 
1.135  mm  (.0447  in)  in  each  plane,  or  with  an  overall 
uncertainty  of  1.605  mm  (.063  in)  in  three-dimensional 
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space.   Temperature  was  measured  at  the  approximate  rate  of 
100  samples  per  second,  well  within  the  constraints  of  the 
thermocouple  time  constant.   Uncertainty  in  the  temperature 
measurements  were  governed  by  the  resolution  of  the  analog- 
to-digital  converter  which  was  .0099  °C.   Uncertainty  in 
the  crossflow  velocity  was  .00625  m/s  (.0205  ft/sec). 
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APPENDIX  B 


MICROCOMPUTER  PROGRAMS 


MAIN 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

MO 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 


rtfiIN_T 

This  program  coordinates  the  entire 
data- talc  i  ng  evolution  for  measuring 
temperature  distributions  in  buoyant 
jets  . 


1 .  Load  all  subprograms 

2.  Input  desired  positions  from  a 
disk  file  of  the  form  :  "RUNEX" 
(value  "XX"  is  the  run  number) 

3.  Move  the  3-D  positioning 
platform  to  each  position 

4.  Align  the  thermocouple  probe  with 
the  jet  streamwise  axis 

5.  Obtain  200  temperatures  at  each 
position  and  compute  the  mean  and 
standard  deviation  (sd) 

6.  Write  to  disk : 

a.  T(C) 

b.  sd 

C.  x  ,  y  ,  z  in  (mm) 

nozzle  centerline  is  (0,0.0) 

and  arm  dimensions  in  inches: 
Shortest  arm  length  -  20.0 
Longest  arm  length   -  24.0 


Probe 
a. 

b. 


Length_probe-4.4375 
Length_arm*20 . 0 
t 

OPTION  BASE  1 

DIM  Coef (12) .X(500).Y(500) ,Z(500) ,Probe_angle(500 ) 

LOADSUB  ALL  FROM  "T_SUBS" 

LOADSUB    ALL   FROM    ••MTR_SUBS•, 

LOADSUB  ALL  FROM  "PROBE_SUBS" 

CALL  Retr ieve_coef (Coef (*) , "motor_coef ") 

2.  Input  desired  positions  from  disk 


1-1 

INPUT 
INPUT 
ASSIGN 


Go 


'Angle  of  arm  relative  to 
'Filename  for  positioning 

'?File4  TO  Filenames 

! 

ENTER  '§iFile4;X(I)  ,Y(I)  ,Z(D 
i 

IF  X(I)O-100  THEN 
IF  X(I)=- 999  THEN 

Probe_angle(I)=Y(I) 

P_angle=Y(I) 
ELSE 

Probe_angle( I) =P_ang le 

1  =  1+1 
END  IF 
t 

GOTO  Go_on 
END  IF 


+  Y-ax  is?" ,Angle_arm 
data7" ,F i lenameS 
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6IC 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

e_ang 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

I  110 

1120 


GOTO  Go_on 

END  IF 
i 

ASSIGN  *File4  TO  » 
Ni  terns*  I- 1 

!3.  Begin  loop  to  take  data  at  each  point 
i 

INPUT  "NAME  OF  FILE  WHERE  DATA  IS  TO  BE  STORED?" ,Fi lenamel $ 

Records  =  <  N i tems»8*7/256 ) +2 

CREATE  BDAT  Fi lenamel $ .Records 

ASSIGN  ©Filel  TO  FilenamelS 

t 

t 

FOR  I_position-1  TO  Nitems 

a.  Move  milling  machine  and  move  the 
thermocouple  probe, 

I-I_pos  i  t ion 
CALL  Move_ldv_to<X(I) , Y< I) , Z< I) .Lengh t_arm ,Length_probe , Angle_arm 
le<I) ,Coef (»)) 


b.  Obtain  mean  temperature  and  sd 
CALL  T_couple<T,,,PRQBE"."C".200,St_dev) 

c.  Measure  ambient  and  jet  temps 


CALL  T_couple<T_ajnbient, "AMBIENT", "C\  10, Sd) 

CALL  T_coup le<  T_noz  z le , "NOZZLE" , "C" , 1 0 , Sd ) 

d.  Write  all  information  to  disk 


OUTPUT  OFilel ;X<I> , Y< I) ,Z( I) 

OUTPUT  @File1 ; T , T_ambien t , T_noz z le 

OUTPUT  i>File1  :St_dev 
NEXT  I_position 
t 

!4.  Close  files 
i 

OUTPUT  ©File! ; -100 

ASSIGN  ©Filel  TO  * 


END 


5.  Terminate  program 

PRINT  "All  done!" 

BEEP 

BEEP 
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PROBE  SUBS 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 


PROBE_SUBS 

This  program  moves  the  temperature  probe 
to  desired  angles  of  deflection, 

NOTE:   Calibration  coefficients  are 
entered  in  SUB  Read_angle, 
beginning  at  line  290. 

CALL  Read_angle(Angle> 

PRINT  "The  probe  is  presently  at";Angle;"  degrees  from  horizontal 
INPUT  "What  is  your  desired  angle  for  the  probe?" , Des l red_ang le 
CALL  Probe_move< Desi red_angle) 

GOTO  120 
END 


SUB  Read_angle( Actual_ang le) 

This  program  reads  the  present  angle 
'of  the  probe  and  returns  it. 

0  degrees  =  horizontal 
'90  degrees  =  vertical,  downward 

A0 ,  A1  and  A2  ■  coefficients  for  a 
second  order  curve  fit  of  mv  vs  angle 
of  deflection  data. 


Coefficients  for  8  June  data  follow 


--23.8657824056449 
=3.09995506283164E-2 
=  -1  .22922305456149E-6 

!  FORMAT  A/D  CARD 
TPUT  723;"CC,1T"  !Clear  Relay  Card 
TPUT  723;"CC,3T"  IClear  A/D  Card 
TPUT  723;"CC,7T,,  !Clear  Digital  Card 
TPUT  723  ;"SF ,3,3,3,1  .25.12T" 
rPUT  723:  "OBJ  ,  10,  IT"  !  CLOSE  RELAY 
TPUT  723;"IP,3T"     !  START  A/D 
TER  72301 ;v 
tual_angle=A0+A1»V+A2»V*V 


A0 
A1 
A2 
» 

QU 

OU 

OU 

OU 

OU 

OU 

EN 

Ac 
SUBEND 
I 
» 
I 
SUB  Probe_move( Desi red_angle) 

This  subprogram  moves  the  probe  to  the 
'desired  angle 


0.  Check  to  see  if  the  angle  is  in 
an  acceptable  range. 

OUTPUT  723:"0P. 1 ,0T"  IClear  Relays 
CALL  Clear_screen 

IF  Desired_angle>90  THEN 
BEEP  3400.1 
BEEP  17  00.1 
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600 

BEEP  3400,1 

610 

PRINT  "Desired  angle  exceeds  90  degrees!!! 

620 

SUBEXIT 

630 

END  IF 

640 

i 

650 

IF  Desired  angle<0  THEN 

660 

BEEP  3400.1 

670 

BEEP  3800.1 

680 

BEEP  3600.1 

690 

PRINT  "Desired  angle  is  negative!!!" 

700 

SUBEXIT 

710 

END  IF 

720 

730 

l.a.  Clear  the  digital  output  card. 

740 

b.  Format  the  A/D  card. 

750 

c.  Close  the  relay  that  corrects 

760 

the  probe  potentiometer  to  the 

770 

A/D  converter. 

780 

790 

OUTPUT  723;,,OP.7,0T" 

800 

OUTPUT  723;"SF,3,3,3,1 .25.3T" 

810 

OUTPUT  723;"CC,1T" 

820 

OUTPUT  723;"0B.1  . 1 0  , 1 T" 

830 

840 

2.  Define  the  acceptable  tolerance 

850 

in  the  angle  (degrees). 

860 

870 

Tolerance- . 5 

880 

890 

3.  Control  loop. 

900 

910  Repea 

t: 

CALL  Read  angle(Actual  angle) 

920 

PRINT  "ANGLE  -"; 

930 

PRINT  USING  "DDD.DD";Actual  angle 

940 

BEEP  Actual_angle-100. .05 

950 

Angle  er ror* < Des i red  angle-Actual  angle) 

960 

i 

970 

IF  ABS(Angle  error )>Tolerance  THEN 

980 

t 

990 

IF  Angle_error>-0  THEN 

1000 

Direction£-"Dou»n" 

1010 

ELSE 

1020 

D  l  rect  ion$-"Up" 

1030 

END  IF 

1040 

1 

1050 

CALL  Motor  go<Di rectionS) 

1060 

» 

1070 

GOTO  Repeat 

1080 

END  IF 

1090 

OUTPUT 

723;"OP,7.0T,, 

1100 

OUTPUT 

723;"CC,1T"   !Clear  relay  card. 

1110 

SUBEND 

1120 

t 

1130 

i 

1140 

I 

1150 

SUB 

Mo 

tor  go(Direc t ion$) 

1160 

IF 

DirectionS="Up"  THEN 

1170 

Lbit=8 

1180 

EN 

D  IF 

1190 

t 
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1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 


IF  Direction$="Doun"  THEN 

Lbit=7 
END  IF 
i 

OUTPUT  723;••aP,7M;2••Lblt;,•T,• 
SUBEMD 


SUB  Clear_screen 

!  Clear  the  CRT. 
t 

OUTPUT  2  USING  "* . B" ;255, 75 

GCLEAfi 
SUBEND 
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3.   MTR  SUBS 


10  !    MTR_SUBS 

20 

30  !   The  following  series  of  subroutines 

40  !  are  utilized  to  calibrate  the  positioning 

60  !  platform  and  ultimately  to  move  the  probe 

70  !  tip  to  desired  positions  within  tne  jet. 

80  ! 

'30  SUB  Draw_flume 

100  ! 

110  !  Draw  the  Buoyant  Jet  Flume  on  the  CRT. 

120  ! 

120  GCLEAR 

1  A0  GRAPHICS  ON 

150  WINDOW  0,48,0,38 

160  LINE  TYPE  1 

170  MOVE  6,6 

180  !  Draw  the  top  view 

190  IDRAW  36.0 

200  IDRAW  0,10 

210  IDRAW  -36,0 

220  IDRAU  0,-10 

230  ! 

240  IMOVE  6,0 

250  IDRAU  0.10 

260  IMOVE  24,0 

270  IDRAU  0.-10 

280  !  label  •'Top". 

290  MOVE  0.10 

300  CSIZE  5 

310  LABEL  "Top" 

320  !  Draw  the  side  view. 

330  MOVE  6,22 

340  IDRAU  36,0 

350  IDRAU  0,13 

360  IDRAU  -36,0 

370  IDRAU  0.-13 

380  IMOVE  6,0 

390  IDRAU  0,13 

400  IMOVE  24,0 

410  IDRAU  0,-13 

420  !  label  "Side" 

430  CSIZE  5 

440  MOVE  0.29 

450  LABEL  "Side" 

460  ! 

470  !  Label  the  picture. 

480  ! 

490  MOVE  1 1 ,35 

500  CSIZE  7 

510  LABEL  "BUOYANT  JET  FLUME" 

520  ! 

530  !  Put  on  the  nozzle. 

540  ! 

550  MOVE  14.22 


560  IDRAU  0,2 

570  IDRAU  .25,0 

530  IDRAU  0.-2 

530  ! 

600  IMOVE  0.-2 

610  CSIZE  3 
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620 

630 

640 

650 

660 

670 

630 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1  100 

1110 

1120 

1130 

1140 

1150 

I  1 60 

1170 

I  180 

1190 

1200 

1210 


LABEL  "nozzle" 

t 

MOVE  14.11 
CSIZE  3 
LABEL  "o" 

Indicate  the  direction  of  flow, 

MOVE  8,11 
IDRAW  2,0 
IDRAW  -.5, -.5 
IDRAW  .5,  .5 
IDRAW  -.5,  .5 
i 

MOVE  8,27 

IDRAW  2,0 

IDRAW  -.5, -.5 

IDRAW  .5,. 5 

IDRAW  -.5.  .5 
i 

SUBEND 


SUB  Calibrate(Filename$> 
t 

OPTION  BASE  1 
DIM  Coef(12) 

Calibrate  the  positioners  on  the 

milling  machine  movement. 

Onto  disk,  write  oat  the  calibration 

coefficients  and  the  hard  boundaries 

that  mast  be  observed! 

This  file  will  be  called  "Motor  coef". 


Position  the  probe  volume  at  the 
wall  of  the  tip  of  the  nozzle.  This 
position  is  (0,0,0).   All  readings 
will  be  in  inches.   Read  all  three 
potentiometers.   Ask  the  user  for  the 
nozzle  outer  diameter  and  compute  the 
zero  position.   Ask  the  user  for  the 
milling  machine  readings. 

Next,  move  the  bed  to  some  new  posi- 
tion using  the  override  switches. 
Take  readings  from  the  pots  and  ask 
for  the  milling  machine  readings. 
Compute  the  calibration  coefficients. 

Move  the  bed  to  each  of  the  extremes 
in  the  X,  Y.  and  Z  directions  using 
the  override  switches  and  have  the 
user  tell  the  computer  when  each  of 
these  boundaries  are  hit.   Enter  each 
of  these  onto  the  disk  file. 
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1220 

1230 

1240 

I  250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

I590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 


D.   Disk  file  "Motor  coef" 


1 . 

x_zero 

x_s lope 

2. 

y_zero 

y_s lope 

3. 

z  zero 

z_s lope 

4. 

x_m i n , 

x  max 

5. 

y_min , 

/_max 

b. 

z_min , 

Z   IDiX 

BEEP 

PRINTER  IS  1 

GCLEAR 

OUTPUT  2  USING  "*  ,B" ;255 . 75 

i 

POTENTIOMETER  CALIBRATION 


PRINT  "I 
PRINT  •' 
PRINT  " 
PRINT  " 
PRINT  •* 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT 
PRINT 
PAUSE 
BEEP  1500 
INPUT  "1 . 
BEEP  2000 
INPUT  "2. 
BEEP  2500 
INPUT  "3. 
BEEP  3000 
INPUT  "4. 


NOTE:  1.  Probe  must  be  horizontal" 

2.  Arm  mast  be  parallel  to 
the  bed  axis 

3.  ALPHA  -  BETA 

•• 

A.  Using  the  override  switches, 
position  the  probe  volume  at 
the  outer  mall  of  the  tip  of 
the  nozzle.  " 

B.  I  will  need  the  nozzle  O.D.  and" 
the  milling  machine  positon. 

(  Hit  <cont>> 


.1 

Nozzle  O.D.  ( inches)?" ,Noz z le_od 
,  .1 

X  (  in)  ,  {+  into  f  lumeP",X_1 
,  •  1 

y  (in),  {+  along  flume  to  the  r lght I?" , Y_1 
,  .  1 

z  (in),  (+  upward}?" , Z_1 


CALL  Read_pot("X",vx_1> 
CALL  Read_pot<"Y",v,y_1  ) 
CALL  Read_pot("Z",Vz_1 ) 
I 
i 

CALL  Clear_screen 
PRINT  "C .  Move  the  milling  machine  to  a  new" 
PRINT  "    position  in  3-D,  by  at  least  5 
PRINT  "    inches  in  each  direction. 

BEEP  3200,  .  1 
INPUT  "X,  Y,  Z  in  inches'?"  ,X_2  ,Y_2  ,  Z_2 

CALL  ReadjDot("X",Vx_2) 

CALL  Read_pot("Y",vy_2) 

CALL  Read_pot("Z",Vz_2> 


C.   Calculate  the  calibration  coefficients 

X_zero--(<Nozzle_od/2)+Vx_l»((X_2-X_1 ) / ( Vx_2-vx_1 ))) 
X_slope=(X_2-X_1 )/(Vx_2-Vx_1 ) 
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1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 


i 

Y_zero— Vy_1-(Y_2-Y_1)/<Vy  2-Vy_1  > 

Y_slope=( Y_2-Y_1 )/(Vy_2-Vy_1 ) 
t 

Z_zero— Vz_WZ_2-Z_1>/<Vz_2-Vz_1  > 

Z_sloPe=(Z_2-Z_1 >/<vz_2-v'z_1 > 


D.  Find  the  physical  boundaries  for  each 
di  rect ion . 

CALL  Clear_screen 
PRINT  "1 .  Move  the  milling  machine  to  the 
PRINT  "    minimum  value  of  'x'.  Hit  <cont>. 
PRINT  "(Away  from  the  flume,  backwards; 
PAUSE 

CALL  Read_pot<"X",v) 

X_min=X_zero+X_s lope»V 
PRINT  "2.  Move  the  milling  machine  to  the 
PRINT  "    maximum  value  of  'x'.   (Towards 
PRINT  "    the  flume).   Hit  <cont>. 
PAUSE 

CALL  Read_pot("X",V) 
X_max=X_zero+X_slope»V 
PRINT  "3.  Move  the  milling  machine  to  the 
PRINT  "    minimum  value  of  'y'.  Hit  <cont> . 
PRINT  "(To  the  left  along  the  flume) 
PAUSE 

CALL  Read_pot("Y'\v> 

Y_min  =  Y_zero+Y_slope*v' 
PRINT  "4.  Move  the  milling  machine  to  the 
PRINT  "    maximum  value  of  'y'.  Hit  <cont>. 
PAUSE 

CALL  Read_pot("Y",V) 

Y_max»Y_zero+Y_slope*V 
PRINT  "5.  Move  the  milling  machine  to  the 
PRINT  "    minimum  value  of  'z'.  Hit  <cont> 
PRINT  "{Downwards}" 
PAUSE 

CALL  Read_pot("Z",V> 
Z_m in=Z_zero+Z_slope*V 
PRINT  "6.  Move  the  milling  machine  to  the 
PRINT  "    maximum  value  of  'z'.  Hit  <cont> 
PAUSE 

CALL   Read^otCZ"  ,V> 

Z_max=Z_zero+Z_s lope»V 


E.  Write  out  the  file  "Motor  coef". 


ON  ERROR  GOTO  Purge_fiie 
Reenter:  CREATE  BDAT  Filenames, 1 
ASSIGN  SFile  TO  Filenames 

OUTPUT  3File:X_zero,X_slope 
OUTPUT  AFile;Y_zero.Y_siope 
OUTPUT  9Fi le:Z_zero,Z_slope 
t 

OUTPUT  *File;X_min,X_max 
OUTPUT  aFile:Y_min, Y_max 
OUTPUT  ■s>File;Z_min.Z_max 
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2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2580 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 


ASSIGN  .SFile  TO  - 
SUBEXIT 
Purge_file:   PURGE  Filename* 
GOTO  Reenter 
» 

SUBEND 


SUB  Read_pot( Directions, Value) 

Read  one  potentiometer  and  return  a  volt- 
age . 

R3  --  Pot  X 
R4  --  Pot  Y 
R5  --  Pot  Z 


IF  Directin$="X"  THEN  Relay-3 
IF  Direction$="Y"  THEN  Relay=4 
IF  Directions-"!"  THEN  Relay-5 


OUTPUT  723:"0P,1 ,0T" 

OUIPUT  723;"0B,1 , " ;Relay ; " , 1 T' 

OUTPUT  723;"IP,3T" 

ENTER  72301 ;Value 
i 


SUBEND 


SUB  Clear_screen 

!   Clear  the  CRT. 
t 

OUTPUT  2  USING  "*,B";255,75 
GCLEAR 
SUBEND 


SUB  11otor(Direction$. Rotations) 

Turn  on  the  motor  in  the  requested  direc- 
tion <x,y,z)  with  the  requestion  rota- 
tion (CM.  CCW). 
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3020 
3030 
3040 
2050 
30G0 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
32S0 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 

3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 


Di  rS=DirectionS 

RotS=Rotat ion$ 


IF 
IF 
IF 
IF 
IF 
IF 
t 


DirS«"X" 
Dlr$-"X" 
Dir«-"Y" 
Dir$-"Y" 
Dir$-"Z" 
Dir$-"Z" 


AND 
AND 
AND 
AND 
AND 
AND 


RotS="CW"  THEN  Lbit-2 
RotS-"CCW"  THEN  Lbit- 1 
RotS-'TU"  THEN  Lbit-4 
Rot$="CCW"  THEN  Lbit-3 
RotS-"CCW"  THEN  Lbit-5 
Rot$«"CW"  THEN  Lbit-6 


OUTPUT  723 ; "OP ,  7 , " ; 2~Lb i t ; "T" 

SUBEND 


SUB  Motor_stop 

Stop  all  motors! 

OUTPUT  723:"OP,7,0TM 

SUBEND 


SUB  Retr ieve_coef (Coef ( ») ,F i lenameS) 

OPTION  BASE  1 

Retrieve  the  potentiometer  calibration 
coefficients  from  a  disk  file  called 
"Motor_coef " .   Place  these  in  an  array 

ASSIGN  SFile  TO  Filenames 
t 

3430    FOR  1-1  TO  12  STEP  2 
ENTER  »File;Coef (I) ,Coef<I+1) 
NEXT  I 
t 

ASSIGN  $File  TO  » 
SUBEND 


SUB  Label_point(X,Y,Z, Symbols) 

Label  a  point  on  the  Flume  diagram 
using  the  symbol  specified. 

PRINT  TABXYd  ,13)  ."(X.Y.Z)  inches  ■ 
PRINT  USING  MDDD.DDD";X,Y.Z 
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3S20 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

)) 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 


IF  Symbol5-"X"  THEN 
t 

MOVE  14+Y.24+Z 

IMOVE  -.2. -.2 

IDRAW  0, .4 

IDRAW  .4,0 

IDRAW  0.-.4 

IDRAW  -.4,0 

i 

MOVE  14+Y.1 l+X 

IMOVE  -.2, -.2 

IDRAU  0,  .4 

IDRAW  .4,0 

IDRAW  0.-.4 

IDRAW  -.4,0 

END  IF 

i 

IF  Symbols-"-*-"  THEN 
i 

MOVE  14+Y.11+X 
IMOVE  0..2 
IDRAW  -.2. -.4 
IDRAW  .4,0 
IDRAW  -.2, .4 

MOVE  14+Y.24+Z 

IMOVE  0.  .2 
IDRAW  -.2, -.4 
IDRAW  .4,0 
IDRAW  -.2,  .4 
END  IF 
SU8END 


SUB  Move_ldv_to<  X , Y , Z ,Length_arm ,Length_probe , Ang le_arm , Ang le_prcbe ,Cc 

DEG 

t 

OUTPUT  723;"CC,1T" 
OUTPUT  723;"SF.3,3,3,1 .25.12T" 
PRINTER  IS  1 

ON  KEY  0  LABEL  "ABORT"  CALL  Stop_all 
OPTION  BASE  1 
Tolerance- . 006 
» 

!   Move  the  probe  to  the  position  indicated 

!   in  (inches)  relative  to  the  nozzle  tip. 
t 

!   Dimensions  are  in  inches  and  degrees 

i 

!   A.1   Load  in  the  calibration  coefficients. 

X_zero-Coef ( 1 ) 
X_slope=Coef <2> 
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4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
43G0 
4370 
4380 
4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

Angle_p 

4670 

Ang le_p 

4680 

Ang le_p 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 


Y_zero=Coef (3) 
Y_slope  =  Coef  <  4) 
Z_zero=Coef (5) 
Z_slope=Coef ( 6) 
X_min=Coef (7) 
X_max:=Coef  (8) 
Y_min-Coef (9) 
Y_max=Coef (10) 
Z_min=Coef (1  1) 
Z_max=Coef (12) 

A. 2  Move  the  probe 

CALL  Probe_moe( Angle_probe) 

B.   Does  (X.Y.Z)  He  within  the  per- 
mitted boundaries^ 

IF  X>X_max  OR  X<X_min  OR  Y>Y_max  OR  Y<Y_min  OR  Z>Z_maj  OR  Z<Z_min  THEN 

BEEP  1700,  .5 
BEEP  2000,  .5 
BEEP  1700,  .5 
BEEP  2000,  .5 

PRINT  "Desired  point  is  oat  of  range!" 
SUBEXIT 
END  IF 

C.   Find  oat  where  the  probe  is  now,  draw 
the  flame  on  the  CRT,  and  label  the 
desi red  posi t ion . 

1.   Sound  warning,  movement  immenent! 

CALL  Clear_screen 

PRINT  TABXYd  ,  12)  ."MOVEMENT  OF  MILLING  MACHINE  IMMINENT!!!" 
FOR  1=1  TO  4 
BEEP  1200,  .1 
BEEP  1700,  .1 
BEEP  2200,  .1 
BEEP  2700,  .1 
NEXT  I 

CALL  Clear_screen 
t 
i 

OUTPUT  723;"CC,1T"    ! CLEAR  RELAY  CARD. 
CALL  Posi t ion( "X" , X_actual , Valu_shaf t ,Length_arm , Length_probe .Ang le_arm . 
robe ,Coef ( *) ) 

CALL  Posi t ion("Y" , Y_actual , Valu_shaf t ,Length_arm ,Length_probe . Angle_arm . 
robe ,Coef (*) ) 

CALL  Pos  1 1  lonC'Z'",  Z_actual  ,  Valu_shaf  t  ,Length_arm  ,  Length _p  robe  .  Ang  le_arm  , 
robe , Coef ( *) ) 
» 
t 

CALL    Draw_flome 

CALL  Label_point (X_actoal , Y_actual , Z_actual ,"+") 
CALL  Label_point(X.Y,Z,"X"') 

D.   Move  each  motor  to  bring  the  error 
between  actual  and  desired  position 
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4770     !       into  tolerance. 

4780     ! 

4790         X_old-X_actual 

4800  X_node:  Xer ror =X-X_actual 

4810         IF  ABS<Xerror)>Tolerance  THEN 

4820  IF  Xerror>0  THEN  Rot$="CCW" 

4830  IF  Xerror<0  THEN  RotS-'TU" 

4840  CALL  MotorC'X"  ,Rot$> 

4850    CALL  Posi t ion<"X" ,X_actual . Valu_shaf t , Length_arm,Length_probe, Angle_ai 

Ang le_probe , Coef ( *) ) 

4860  "  CALL  P Lot_path< X_old , Y_actual , Z_actual , X_actual , Y_actual . Z_a( 

al) 

4870  X_old-X_actual 

4880  GOTO  X_node 

4890         END  IF 

4900  ! 

4910  ! 

4920         Y_old-Y_actual 

4930  Y_node :  Yerror= Y-Y_actual 

4940         IF  ABS(Yerror)>Tolerance  THEN 

4950  IF  Yerror>0  THEN  Rot$-"CCW" 

4960  IF  Yerror<0  THEN  Rot$="CU" 

4970  CALL  Motor ( "Y" ,Rot$> 

4980    CALL  Posi tionC'Y" , Y_actual , Valu_shaf t ,Length_arm ,Length_probe , Angie_a 

Angle_probe,Coef<">> 

4990  "  CALL  Plo t_path(X_actual , Y_old , Z_actual ,X_actoal , Y_actual , Zi 

al) 

5000  Y_old-Y_actual 

5010  GOTO  Y_node 

5020         END  IF 

5030  ! 

5040  ! 

5050         Z_old=Z_actual 

5060  Z_node :  Zerror  =  Z-Z_actuLal 

5070         IF  ABS(Zerror)>Tolerance  THEN 

5080  IF  Zerror>0  THEN  Rot$-"CH" 

5090  IF  Zerror<0  THEN  Rot$-"CCW" 

5100  CALL  Motor("Z".Rot$) 

5110    CALL  Posi t ion("Z" , Z_actual , Val u_shaf t ,Length_arm ,Length_probe . Angle_a 

Ang le_probe , Coef ( ») ) 

5120  "  CALL  Plot_path(X_actual ,Y_actual . Z_old , X_actual . Y_actual , Z_a 

al) 

5130  Z_old-Z_actaal 

5140  GOTO  Z_node 

5150         END  IF 

5160  ! 

5170  ! 

5180    CALL  Motor_stop 

5190  ! 

5200  ! 

5210    FOR  1-1  TO  4 

5220       BEEP  2400,  .2 

5230       BEEP  4800,  .2 

5240    NEXT  I 

5250  ! 

5260  SBEND 

5270  ! 

5280  ! 

5290  ! 

5300  ! 


5310  ! 

5320  ! 

5330  SUB  Posi tion( Direct ion$, Value, Valu_shaf t ,Length_arm  ,Length_probe  .  Ang  le_3.rm 

Ang le_probe , Coef ( * ) ) 

5340    ! 

OPTION  BASE  I 

DEG 

t 

!   Return  the  position  (inches)  for  the 
!     appropriate  direction  relative  to  the 
!     noz  z le  tip. 


5350 
5360 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
54G0 
5470 
5480 
5490 
5500 
5510 
5520 
5530 
5540 
5550 
5560 
5570 
5580 
5590 
5600 
5610 
Ang le_arm) 


CALL  Read_pot (Directions, Voltage) 

X_zero  =  Coef  ( 1 ) 
X_slope=Coef (2) 
Y_zero=Coef (3) 
Y_slope=Coef (4) 
Z_zero=Coef (5) 
Z_s lope=Coef (6) 


Dir$-Direction$ 

IF  DirS«"X"  THEN 

Valu_shaf  t-X_zero+X_s lope* Vol tage 

Value- Valu_shaf t-Length_arm*COS( Angle_arm)»TAN( 45-Angle_arm/2 . 0 ) 

END  IF 
t 

IF  DirS^'Y"  THEN 

Valu_shaf  t=Y_zero+Y_s lope»Vo I tage 

Value-Valu_shaf t+Length_probe*( 1 . -C0S( Angle_probe) )+Length_arm*C0S( 


5620 
5630 
5640 
5650 
5660 
5670 
5680 
5690 
5700 
5710 
5720 
5730 
5740 
5750 
5760 
5770 
5780 
5790 
5800 
5810 
5820 
5830 
5840 
5850 
5860 
5870 
5880 


END  IF 

I 

IF  DirS-"Z"  THEN 

Valu_shaf  t=Z_zero+Z_slope*Vol tage 
Value=*Valu_shaf  t-Length_probe*SIN(Angle_probe) 

END  IF 
t 

SUBEND 


SUB  Stop_all 

!   STOP  ALL  MOTORS  AND  QUIT 
» 

GCLEAR 

CALL  Motor_stop 

PRINT  "MOTOR  CONTROL  ABORTED!!!" 
PRINT  "(HIT  <C0NT>  TO  CONTINUE)" 
PAUSE 
SUBEND 
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5390 
5900 
5910 
5920 
5930 
5940 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
6060 


SUB  Plot_path(X1 ,Y1 ,Z1 .X2.Y2.Z2) 

Plot  the  path  of  the  probe  on  the  flume 
diagram  as  the  motors  move  the  bed. 

Lower  plot  followed  by  upper  plot. 


PRINT  TABXYd  ,  13)  ,"<X,Y,Z> 

PRINT  USING  "DD.DDD";X2,Y2 
MOVE  14  +  Y1 ,24+Z1 
DRAW  14+Y2.24+Z2 

MOVE  14+Y1 ,11+XI 
DRAW  14+Y2.11+X2 
t 

SUBEND 


inches 
Z2 
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4.   T  SUBS 
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!    T_SUBS 

SUB  T_coupIe( Temperature ,Choice5,ScaleS,No_readings,Stdev) 

SUBPROGRAM  T_C0UPLE 

by  Bill  Culbreth 
19  April  1984 

PURPOSE:   This  program  is  designed  to 
read  type  T  or  E  thermocouples  and 
return  the  actual  temperature  to 
the  calling  routine. 

Temperature  ~  Temperature  from  the 

thermocouple  in  degrees  F  or  C  from 

the  thermocouple  identified  by 

"Cho  ice$"  . 
ChoiceS  --  Thermocouple  choice,  current- 
ly:  "AMBIENT",  •"NOZZLE",  or,  "PROBE". 
ScaleS   --  "F"  for  Fahrenheit  or  "C"  for 

Celsius,  or  "H"  for  histogram  in  "C" . 
No_readings  ~  Hou  many  readings  of  the 

same  thermocouple  should  the  routine 

take'' 
Stdev  --  The  standard  deviation  of  the 

temperature  for  the  indicated  number 

of  readings  in  the  units  given  by 

"ScaleS" . 


1.  Open  all  relays  and  initialize  the 
A/D  conver  ter . 

OUTPUT  723 :"CC, IT- 
OUTPUT  723;"SF.3.3,3,1 .25.12T" 
» 

OUTPUT  723;"0P,1 ,0T" 

2.  Close  the  chosen  relays. 

a.  Ambient  T  ~  Type  T,  relays  6,8. 

b.  Nozzle   "  --  Type  T,  relays  7,8. 

c.  Prooe    T  --  Type  E,  relay   9. 

IF  Scale$="H"  THEN 

!  Plot  a  histogram  using  Celsius. 

HistogramS="YES" 

ScaleS="C" 
ELSE 

HistogramS="NO" 
END    IF 
I 
i 

IF   ChoiceS-"AMBIENT"    THEN 
TypeS="T" 

OUTPUT    723;"0B,1 .6,1 ,8. IT" 
END    IF 
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IF  ChoiceS-"NOZZLE"  THEN 

Type$="T" 

OUTPUT  723  :'"QB,  1  ,7.1  ,3,  IT" 
END  IF 
t 

IF   ChoiceS-'TROBE"    THEN 

Type$="E" 

OUTPUT    723;"OB,1 .9, IT" 
END    IF 
i 

3.   Take  an  A/D  conversion  and  convert 
into  temperature. 

Sum=*0 
SumJ-0 

WAIT  1 

IF  Histogram$-"YES"  THEN  GOSUB  Set_up_histo 
) 

! 

FOR  I« I  TO  No_readings 
OUTPUT  723;,,IP,3T" 
ENTER  72301 :A 

IF  Type$-"T"  THEN 

A-A/1000 
END  IF 
IF  Type$-"E"  THEN 

A-A/1000 
END  IF 
!  PRINT  "Vdnv')  -  ";A 

GOSUB   Convert_t 
!  BEEP   A-100, .01 

Sum-Sum+A 

Suml -Sum1 +A-A 
i 

IF  Histogra/n$«MYES"  THEN  GOSUB  Plot_point 
i 

NEXT  I 

Tempe ra t ur e- Sum/ No_ readings 

( 

i 

IF  No_readings-1  THEN 

Stdev=0 
ELSE 

Stdev-SQR(ABS((Sam1-No_readings*Temperature'2)/(No_readings-l > > > 

END  IF 

4.   Open  all  relays. 

OUTPUT  723 ;"0P , 1 ,0T" 
t 

SUBEXIT 


Convert  t:   ! 


!  Thus  subroutine  converts  (mV) 
!  from  a  thermocouple  into  Temp- 
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!  erature. 


IF  Type$-"T"  THEN 

A-2.5661297E+1*A-6.1954869E-1»A»A+2.2181644E-2»A-3-3.55009E-4»A~4 

END  IF 

t 

IF  Type$-"E"  THEN 

A-1 .7022525E+1*A-2.2097240E-1*A-A+5.4809314E-3»A'3-5.7669892E-5»A~4 

END  IF 

Fix  the  scale. 

IF  Scale$-"F"  THEN  A-1.8-A+32 
RETURN 


Set_up_histo : 
i 


Set  up  a  histogram  of  temperature 
versus  number  of  counts. 

1.  Zero  out  the  Height(*)  array. 

DIM  Height(203) 
i 

FOR  1-1  TO  202 

Height<I)-0 
NEXT  I 
i 
» 

DUMP  DEVICE  IS  701 

GINIT 

OUTPUT  2  USING  "* ,B" ;255 , 75 

GRAPHICS  ON 

FRAME 

WINDOW  -100,100,-10,100 

MOVE  -65,92 

pry  ~7C      7 

LABEL  "TEMPERATURE  HISTOGRAM" 
AXES  25,10,0,0,4,5,3 
PEN  -1 
MOVE  0,-10 
DRAW  0,0 
MOVE  0,90 
DRAW  0.100 
PEN  1 
\ 

!   Take  10  temperature  readings  to  get 
!     the  scale. 


Sum-0 
Sum1=0 
FOR    1-1    TO    10 

OUTPUT    723;"IP.3T" 
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ENTER  72301  :A 
NEXT  I 
t 

FOR  1=1  TO  100 

OUTPUT  723;•'IP,3T,• 

ENTER  72301  :A 

IF  TypeS="T"  THEN 
A-A/1000 

END  IF 

IF  Type$-"E"  THEN 
A=A/1000 

END  IF 

GOSUB  Convert_t 

S  am  *  S  am  +  A 

Saml =  Sum1 +A»A 
NEXT    I 

T_mean=Sum/1 00 

Sd-SQR(ABS(Sam1-100»T   mean"2)/99> 


Change  the  window  to  extent  from 
-4»Sd  to  +4»Sd. 

MOVE  -12.-7 

CSIZE  4 

LABEL  USING  "DDD .DD" ; T_mean 

MOVE  15.-7 

LABEL  "C" 

» 

MOVE  -87,-7 

LABEL  USING  "DDD . DD" ; T_mean-3*Sd 

MOVE  63,-7 

LABEL  USING  "DDD . DD" ; T_mean+3-Sd 

i 

» 

MOVE  5.47 

LABEL  "50" 
t 

» 

!   3.  Calculate  the  window  temper- 

!      ature  interval . 
» 

interval-4»Sd/100 
t 

T_min-T_mean-4»Sd 
T_max=T_mean+4»Sd 
» 

WINDOW  T_min,T_majc,-10,  100 


RETURN 


Plot_point:   ! 


Plot  each  temperature  point  as 
received  on  the  histogram. 

Freq=((A-T  min)/(T  max-T  m  in ) )»3000  +  1 000 
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2410 

BEEP  Freq 

,  .01 

2420 

J  =  0 

2430 

f 

2440 

FOR  T-T  nun  TO  T  max  STEP  Interval 

2450 

J-J+1 

2460 

IF  fl<  = 

T  THEN 

2470 

Height(J)=Height(J)+1 

2480 

GOTO  Continue 

2490 

END  IF 

2500 

NEXT  T 

2510 

2520 

?   Draw  the  point. 

2530 

2540 

Cont inue : 

t 

2550 

t 

2560 

» 

2570 

MOVE  P 

.Height(J) 

2580 

CSIZE 

4 

2590 

LABEL 

"  .  " 

2600 

2610 

!   Label 

the  number  of  points  and 

2620 

!     the 

temperature . 

2630 

2640 

PRINT 

TABXY(2,3>;"T(C)-";A 

2650 

PRINT 

TABXY<2.5> ; "Sample  *" ;I 

2660 

f 

2670 

RETURN 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

SUBEND 
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T_CAL 

This  program  records  the  output  from  the 
ambient,  probe  and  nozzie  themocouples 
for  calibration  purposes.   The  data  is 
printed  as  well  as  transmitted  to  a  data 
file. 

Data  is  recorded  in  the  following  order: 

a.  Ambient  temperature 

b.  Probe  temperature 

c.  Nozzle  temperature 

Place  the  ambient,  probe  and  nozzle 
thermocouples  into  a  bucket  of  warm  water 
execute  this  program,  then  add  ice. 


L0ADSU8  ALL  FROM  "T_SUBS" 

Identify  the  BDAT  file 
NPUT  "NAME  OF  FILE  WHERE  DATA  IS  TO  BE  STORED?" ,  Fi  lenamel  S 

Identify  the  number  of  data  points  desired 

NPUT  "NUMBER  OF  DATA  POINTS  DESIRED?" , N i terns 

Identify  the  number  of  samples  per 
thermocouple  per  data  point 

INPUT  "NUMBER  OF  SAMPLES  PER  DATA  POINT?" .Sample 

Initialize  a  counter 

I 

Create  a  data  file 

Records- <Nitems»8»3/256>+2 

CREATE  BDAT  F i lenamel S .Records 
ASSIGN  9Filel  TO  FilenamelS 

Take  data  and  print  results 

CALL  T_coupie(T_ambient,"AMBIENT","C",Sample,Sd> 

BEEP  T*50.  .05 

PRINT  "TAMB,SD=":T_ambient,Sd 

CALL  T_coup le< T , "PROBE" , "C" .Sample ,St_dev> 

BEEP  T-50,  .05 

PRINT  "TPROBE.SD=";T,St_dev 

CALL  T_couple(T_nozzle, "NOZZLE", "C", Samp le.Sd) 

BEEP  T-50,  .05 

PRINT  "TNOZL.SD=";T_nozzle,Sd 

t 

*   Send  data  to  the  data  file 
t 

OUTPUT  3Filel:T  ambient, T.T  nozzle 

!  ~  - 

!   Test  to  see  if  finished 
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IF  I=Ni terns  THEN 
GOTO  760 

ELSE 

r-i+i 

GOTO  460 

END  IF 

Close  the  data  file 

OUTPUT  SFilel : - T  00 
ASSIGN  SFilel  TO  * 

Alert  the  user. ..the  job  is  completed 

PRINT  "ALL  DONE" 

3EE  = 
3EE? 
EM 
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PROBE  CAL 
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PR0BE_CAL 

This  proqram  will  provide  the  data 
necessary  to  calibrate  probe  deflection 
if  used  in  the  following  manner: 

a.  Switch  OFF  probe  actuation  power 
at  the  probe  base. 

b.  Swing  the  probe  arm  out  of  the 
tank  such  that  the  probe  is 
accessab le . 

c.  Attach  the  calibration  panel  to  the 
probe  assembly,  taking  care  not  to 
damage  the  glass  probe. 

d.  Run  the  program... it  will  ask  for 
the  desired  position  in  mV .   The 
following  guidelines  apply: 

i.  If  it  is  desired  to  lower 
the  probe,  type  the  extreme 
value  9000  (anything  >4600 
will  work ,  but  9000  is  a 
quick  and  easy  number  to 
enter ) . 

11.  If  it  is  desired  to  raise 

the  probe,  enter  the  extreme 
value  90  (anything  less  than 
940  wi 1 1  work ) . 

e.  The  program  will  next  ask  which  bit 
is  selected  to  be  "high".   The 
following  guidelines  apply: 

l.  If  it  is  desired  to  lower  the 
probe ,  enter  7 . 

ii.  If  it  is  desired  to  raise  the 
probe ,  enter  8 . 

f.  Switch  ON  probe  activation  power  and 
when  the  probe  reaches  a  desired' 
degree  of  deflection,  switch  the 
activation  power  OFF  and  record  the 
mV  value  printed  on  the  screen. 

g.  Repeat  steps  (d)  through  (f)  until 
sufficient  data  is  collected.   Load 
this  data  into  the  program  "POLYFII" 
and  request  a  second  order  fit 

( let  X=angle  and  Y=mV> . 
Enter  the  coefficients  derived  by 
this  program  into  the  appropriate 
location  in  the  program  "PROBE_SUBS" 


MOTE 


The  relays  will  be  activated  on 
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digital  low!   When  the  machine  boots 
up  (hp-9826),  all  relay  lines  are 
high  (+5V).   The  instructions  belou 
will  drop  the  voltage  to  zero. 

INPUT  "What  is  the  desired  position?  (mV )"  ,  Vol  tage 
PRINT  "Desired  position  (mv1)  =";Voltage 

Read  the  actual  motor  position. 

If  the  desired  position  is  BELOW  the 
actual  position,  then  tell  the  motor 
to  move  UP  . 
If  it  is  ABOVE,  then  tell  it  to  go  DOWN. 


INPUT  "WHICH  BIT  DO  YOU  WANT  HIGH?" ,Lbit 
» 

OUTPUT  723;"OP,7,0T"  !CLEAR  ALL  D/0 

OUTPUT  723;"0P,7,";2"Lbit;"T" 

t 

OUTPUT  723;"CC,1T"  !CLEAR  A/D  CARD 
» 

OUTPUT  723;"SF,3,3,3,1 .25.12T"  !F0RMAT  A/D 
OUTPUT  723;"0B,1  ,10,1T"  ICLOSE  THE  RELAY 

!THAT  CONNECTS  THE 
!A/D  CARD  TO  THE 
POTENTIOMETER 
OUTPUT  723:"IP,3T"  !START  A/D  CONVERSION 
ENTER  72301 ;A      'ENTER  A/D  VALUE  INTO  A. 
DISP  A 
BEEP  ABS(A), .01 


The  following  IF  stops  the  motor  when 
it  reaches  it's  limits. 

IF  (Lbit-7  AND  A>4600)  OR  (Lbit-8  AND  A<940)  THEN 
OUTPUT  723;"OP,7,0T" 
BEEP 
PRINT  "ALL  DONE!!'" 

END  IF 
GOTO  890 
END 
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7.   MOTOR  CAL 
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M0T0R_CAL 

PURPOSE:  Calibrate  the  potentiometers 
used  with  the  motors  on  the  milling 
mach ine . 


OPTION  BASE  1 

DIlj  Coef  (12) 
t 

OUTPUT  723;,'SF.3.3.3.1  .25. 12T" 

LOADSUB  ALL  FROM  "MTR_SUBS" 
LOADSUB  ALL  FROM  "PROBEJSUBS" 
F i leS="motor_coef " 
t 

CALL  Calibrate(File$> 

CALL   Retneve_coef  (Coef  (*)  ,Fi  leS) 

BEEP   2400,1 
i 
» 

PRINT  "Uhat  (x.y.z)  position  do  you  wish  to" 
PRINT  "   move  to"?  (inches  relative  to  noz-" 
PRINT  "   zle)" 

INPUT  X.Y.Z 
» 

CALL   Move_ldv_to(X . Y , Z ,Length_arm ,Length_probe , Angle_arm ,Probe_angle ,Coe 

t 

BEEP 

GOTO   280 
END 
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lGAD_XYZ 

by  Bill  Culbreth 

30  April  1984 


PURPOSE:  'This  program  will  allow  the 
user  to  enter  desired  (x.y.z)  posi- 
tions of  the  milling  machine  relative 
to  the  nozzle.   The  values  will  be 
stored  on  disk  to  be  utilized  later 
by  MAIN_T. 

HM  X<500) ,Y(500) ,Z(500) 

1.  Input  the  file  name. 

G03UB  Clear_screen 

PRINTER  IS  1 

PRINT  "Input  the  file  name." 

PRINT  "{  I  suggest  'RUNXX'  where  'XX'" 

PRINT  "   is  the  run  number.  5" 

PRINT 

INPUT  Filenames 
i 

2.  Begin  inputting  data. 


GOSUB  Clear_screen 
t 

PRINT  "Do  you  wish  to  append  a  previous  data  file?' 

INPUT  Answeri 

I 

IF  Answer$="YES"  THEN 

INPUT  "Previous  file  name?" , 01d_f l le$ 
ASSIGN  SFilel  TO  01d_fileS 

t 

1  =  0 
t 

'i-I+l 

ENTER  SFilel ;X(I) ,Y(I) ,Z(I) 

OUTPUT  701;"I,X,Y,Z=";I,X<I).Y(I).Z(I) 

IF  X(IX>-100  THEN  GOTO  Loopl 

ASSIGN  SFilel  TO  » 

PURGE  01d_file$ 

Count = I- 1 
ELSE 

Count="0 
END  IF 

t 

GOSUB  Clear_screen 

BEEP 


Loopl 


PRINT  "1.  Input  the  desired  position  in 

PRINT  "  inches  as  X.  Y,  and  Z." 

PRINT 

PRINT  "2.  Terminate  input  by  entering   -100 
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PRINT  "2, 
PRINT  " 
PRINT  "3 
PRINT  " 
PRINT 


Terminate  input  by  entering  '  -  1  0 0  ' ' 

for  X,  Y.  and  I." 

If  you  wish  to  enter  an  orientation  angle," 

enter  ' -999 , or len tat  ion , 0 '  . " 


Begin:  ! 


Count*Count+1 

PRINT  "Item  *":Count 

INPUT  "(X.Y.Z)  in  inches'?"  ,  X  (Count ),  Y(Count)  .  ZCCount ) 

OUTPUT  701  ;"  I ,  X,  Y  ,  Z-"  ;  Court  t,X<  Count)  .Y  (Count)  ,Z(Count) 

IF  X(Count)<>-100  THEN  GOTO  Begin 

End_count=Count-1 
t 

All  data  points  have  been  entered. 

a.  Set  up  new  softkeys. 

b.  Explain  softkeys. 

GOSUB  Clear_screen 
i 

PRINT  "SOFTKEY  LABELS:" 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
» 


0  --  WRITE  data  to  disk ." 

2  --  EDIT  out  bad  data." 

4  --  HARD  copy  the  data  on  printer 

6  --  LIST  data  on  the  CRT." 

8  --  STOP  terminates  the  program." 


ON  KEY  0  LABEL  "WRITE"  GOSUB  Wnte_data 
ON  KEY  2  LABEL  "EDIT"  GOSUB  Edit_data 
ON  KEY  4  LABEL  "HARD"  GOSUB  Hard_copy 
ON  KEY  6  LABEL  "LIST"  GOSUB  List_data 
ON  KEY  8  LABEL  "STOP"  GOTO  Terminate 
GOTO  900 


Clear_screen :    ! 
t 

!   Clear  the  CRT  display. 
t 

OUTPUT  2  USING  "*  ,B" ;255 .75 

GCLEAR 
RETURN 


Hard_copy:  \ 
i 


!  Print  out  all  data  to  the  printer 
i 

PRINTER  IS  701 
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1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 


PRINT  "   Count 

PRINT  "   

PRINT 


X(  incnes ) 


y(count>      Z(count)' 


RETURN 


FOR  1=1  TO  End_count 

PRINT  I.X(I) ,Y(I),Z(I) 
NEXT  I 
« 

PRINTER  IS  1 


list_data:  ! 

List  data  to  the  CRT. 

GOSUB  Clear_screen 

PRINT  "There  are" ; End_count ;  "  data  points." 

PRINT 

INPUT  "Which  point  do  I  start  with?", Start 

INPUT  "which  point  do  I  end  u i th?" ,Ena_aata 
» 

FOR  I-Start  TO  End_data 

PRINT  "I.X.Y.Z-" ;I,X(I),Y(I) ,Z<I> 
NEXT  I 
RETURN 


Edit_data:   ! 

Edit  data. 

GOSUB  Clear_screen 
I 

INPUT  "Uhich  data  point  do  you  wish  to  alter7", I 

i 

PRINT 

PRINT  "For  point  *";!;",    (X.Y.Z)  where:" 

PRINT  X(I)  ;",■•  ;YXI)  ;",";Z(I> 

INPUT  "Type  in  the  new  values : " , X( I) , Y( I ) , Z< I) 
RETURN 


Urite_data:   ! 

f 


Write  data  out  to  a  file  on  disk 
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1810 
I320 
1830 
1840 
1850 
1860 
1370 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
I970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 


X(End_count+1 > -- 1 00 

Y<End_count+1 )--100 

Z<End_count+1 > -- 1 00 
t 

Max_data-3»(End_coant+1 ) 
t 

File_size-INT(Max_data«8/256>+1 
t 

CREATE  BOAT  Fi lenameS.Fi le_si ze 

ASSIGN  «File  TO  Filenames 
t 

FOR  1-1  TO  End_count+1 

OUTPUT  @File;X<I),Y<I).Z(I) 
NEXT  I 
i 

GOSUB  Clear_screen 
BEEP  2400. .3 

PRINT  "File  "-.Filenames;"  has  been  stored!' 
ASSIGN  SiFile  TO  » 
RETURN 


Terminate:  ! 

GOSUB  Clear_screen 

PRINT  "NORMAL  TERMINATION  OF  PROGRAM!" 
END 
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SEND  DATA 


10 
20 
30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

ISO 

160 

170 

!80 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570< 

530 

590 

600 


SEND_DATfl 

To  vAX.IBM.TRS-80. 

HP-9826  TERMINAL  PROGRAM 
[REGUIRES  BINARY  ENHANCEMENT  PROGRAM 
"BEB"!] 

JUNE  30.  1982 
updated  1/5/83 
updated  1/16/84 

BILL  CULBRETH 


Sc-9    !  RS-232  IS  SELECT  CODE  9. 
PRINTER  IS  1   !  PRINTER  IS  CRT. 
Pr-1  !  DEFAULT  PRINTER  IS  CRT 

Prmter_choice-701  !  MY  PRINTER  IS  701 
Bits-7        !  BITS  PER  CHARACTER 
Duplex-0       !  FULL  DUPLEX 
Baud-300      !  BAUD  RATE 
Computer-1     !  ASSUME  IBM  COMPUTER 


(300  BAUD,  IBM  assumed." 
Load  the  binary  program  BEB  first" 
unless  you  have  BASIC  2.0" 
SET  MODEM  ON  <FULL  DUPLEX>  )" 


OUTPUT  Pr 
OUTPUT  Pr 
OUTPUT  Pr 
OUTPUT  Pr 
OUTPUT  Pr 
» 

DIM  NameS [ 200 ],Hp_file$[ 301 ,Aa<l500) .NumbS(30] 

INTEGER  Isend 

i 

CONTROL  Sc,3;Baud 

CONTROL   Sc,4;Bits-5+4    !    BITS/CHAR   A    'STOP   BITS 

I 

» 

fo_dislc-0 

Datadump=0 

I_data=1 

1  =  1 

J-1 

K-1 

L-1 

ON  ERROR  GOTO  Errors 


ON  KEY  0  LABEL 
ON  KEY  5  LABEL 
ON  KEY  6  LABEL 
ON  KEY  7  LABEL 
ON  KEY  8  LABEL 


"Line  Mode"  GOTO  Line_mode 
"Terminal"  GOTO  Terminal 
"To  Crt"  GOTO  Pr_crt 
"To  Prt"  GOTO  Pr_prt 
"DATA"  GOTO  Data_dump 


Line_mode:     ! 

OUTPUT  Pr;"(LINE  RECEPTION  MODE}" 
Begin:  STATUS  Sc  .  1 0 ; Y  !  CHECK  FOR  FULL  BUFFER 

IF  BIT(Y,0)-0  THEN  GOTO  Begin 
t 

!   RECEIVE  ROUTINE. 
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S10    ! 

S20  Receive: 

B  =  A 


STATUS  Sc,6;A 


S30 

640 

650 

660 

670 

680    ! 

690   ! 

700    ! 

710- 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950  . 

960  Data_dump: 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110  GOTO  Begin 

1120 

1130 

1140 

1150 

1160  Errors 

1170 

1  180 

1130 

1200 


OUTPUT  Pr  USING  "*  , A" : CHR$<B> 

IF  B-63  AND  Datadiunp«1  THEN  GOTO  Data_dump 

IF  B-13  AND  Computer-3  THEN  OUTPUT  Pr:CHR$(13) 

GOTO  Begin 

TRANSMIT  ROUTINE. 


IF  Duplex-0  THEN 

IF  NUM<KeyS)<>255  THEN  OUTPUT  Pr  USING  "#,A";Key$ 

IF  NUM<Key$>-255  THEN  OUTPUT  Pr;"  " 
END  IF 

IF  Computer-1  AND  NUM<Key$>-8  THEN  Key$-CHR$(64) 
t 

!  the  previous  line  gives  an  * 
!  for  a  backspace  for  the  IBM. 
t 

IF  Computer-5  AND  NUM(Key$)-8  THEN  Key$-CHR$( 1 27) 
i 

!   THE  VAX/VMS  REQUIRES  A  DELETE 

!   SYMBOL  FOR  A  BACKSPACE. 

t 

IF  NUM(Key$>-255  THEN  Key$-CHR$< 13) 
OUTPUT  Sc  USING  "# ,A";Key$ 
GOTO  Begin 


DATA  FILE  OUT  TO  THE  HOST  COMPUTER. 


IF  I_data-1  THEN  GOSUB  0pen_file 
I 

IF  Datadump-0  THEN  GOTO  Begin 
IF  Computer-1  THEN  WAIT  .3 
?  wait  for  the  slow  IBM. 
BEEP  1000+RND-1500 , .05 
OUTPUT  P^•,,•A<••;I_data;••)-',; 
OUTPUT  Pr ;Aa(I_data) 
GOSUB  Send_nujnber 

IF  Aa(I_data)— 200  THEN 
I_data=1 
Datadmnp-0 

END  IF 

I  data-I  data+1 


ERROR  HANDLING  SUBROUTINE 


OFF  ERROR 
Close_f  ile— 200 

!   FIRST.  END  OF  FILE  ERROR. 
IF  ERRN-59  THEN 
Aa(I>— 200 
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1210 

1220 

1230 

1240 

1250 

12S0 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

13S0 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1550 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 


GOTO  2000  !  RETURN  AFTER  ERROR. 
END  IF 
i 

IF   ERRN059    THEN   OUTPUT   Pr:"<error    #";ERRN;"    generated .  >" 

IF   ERRN=54    THEN   OUTPUT   Pr;"(FILE    <" ;Hp_f i US ; ">    ALREADY    THERE)" 

IF  ERRN-54  THEN  GOTO  Created 

IF  ERRN=56  THEN  OUTPUT  Pr;"<FILE  <" ;Hp_f i leS ; "  IS  NOT  ON  DISK.)' 

ASSIGN  SFile  TO  * 
GOTO  Line_mode 
t 


OUTPUT  TO  CRT, 


Pr_crt:   Pr-1 

GOTO  Line  mode 


!   OUTPUT  TO  PRINTER, 
i 

Pr_prt:   Pr-Pr inter_choice 

GOTO  Line_mode 

» 

i 

j 

Terminal 


CHANGE  THE  TERMINAL  CHARACTERISTICS, 


Baud  Rate  -";Baud 
Bits/Char  =":Bits 
Duplex    -";Duplex 
[  1-fuH.O-half  ]" 
Computer   -";Computer 
CIBM-1 ,  VAX/UNIX-2 
TRS-80-3,  Cyber-4 


GO 
i 
t 
I 

Op 


OUTPUT  Pr;"     1 , 

OUTPUT  Pr;"     2, 

OUTPUT  Pr;"     3, 

OUTPUT  Pr;"' 

OUTPUT  Pr;"     4 

OUTPUT  Pr;" 

OUTPUT  Pr;" 

OUTPUT  Pr;"  " 

INPUT  "Change  which  one?", Which 

IF  Which- 1  THEN  INPUT  "To?", Baud 

IF  Which=2  THEN  INPUT  "To?",8its 

IF   Which-3    THEN    INPUT   "To?", Duplex 

IF  Which-4  THEN  INPUT  "To?" .Computer 

IF  Computer-1  THEN  Duplex-0 

THEN   Duplex=0 

THEN   Bits-8 

THEN   Duplex-1 


vax/vms-5 ]' 


IF  Computer-3 
IF  Computer*^ 
IF  Computer-5 
TO  Line  mode 


en_file:   ! 

!   Open  a  file  to  read  data  from 

!      disk. 

t 

Datadwnp=1 
t 

INPUT  "Is  this  LDV  data?  { 1  -  YES) "  ,Ldv$ 
IF  Ldv$="1"  THEN 

INPUT  "Experiment  #?" .Exper imentS 
ELSE 

OUTPUT  Pr;"Data  file  out  of  HP  to  host." 

INPUT  "File  name?", Hp_f i lei 
END  IF 
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1810 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950- 

1960 

1970 

1980 

1990« 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 


1820     IF  LdyS-'T"  THEN 
Hp_f  ile$  =  ExpenmentSa"_RESULT" 
END  IF 

Read  the  file  off  of  disk . 

ASSIGN  SFile  TO  Hp_file$ 

1-1 

Check-0 

BEEP 

BEEP 

OUTPUT  Pr;M {Working  on  file  <" ;Hp_f i le$;"> . }' 


ENTER  @File;Aa<I> 
Check-Aa(I) 

1  =  1  +  1 

» 

ASSIGN  OFile  TO  • 
Datadmnp-  1 
RETURN 
i 
i 
j 
Send_niunber :  ! 

!  SEND  A  NUMBER  ONE  CHARACTER  AT 
!  A  TIME  TO  THE  HOST  COMPUTER. 
f 

Nujnb$  =  VAL$<Aa<  I_data)  > 

Length=LEN<Numb$) 

i 

IF  (<Ldv$-"1")  AND  <I_data>13>>  THEN 

Posit  =  POS<Ntunb$,"."> 

IF  (PositOO)  THEN  Length-Posit+2 
END  IF 
I 

FOR  1-1  TO  Length 

Numeric«NUM(Numb$tI,I]> 

OUTPUT  Sc  USING  "* ,AM ;Numb$[ I , I] 
NEXT  I 
t 

OUTPUT  Sc  USING  "*  ,  A"  ;CHR$( 13) 
RETURN 
i 
i 
i 
END 


98 


APPENDIX    C 
MAINFRAME    PROGRAMS 


1 .       TCAL 


c 

c 

C        TCAL 

C 

C  PURPOSE:         THIS     PPCGwAM     PLaCcS     THERMOCOUPLE     CALIBRATION     data 

C  SfcCEIVFO     FROM     THE     HP-<3fl26     "  I C  &OCOMPU  T£P      INTO    A 

C  "ORE      <OHKAeLE     FORMAT     FOR     USE      IN      THE     PROGRAM      TF I T . 

C 

c 

O I  MENS I  ON  TPROBtt  2001  .TAMe(200l  .  TNOZZ  <  2  00  1 
I  =  I 
10   REAQI  07 ,* I   TAMR(  I  | 

[F(  r AMH (  I  |  .EO.-200. 0  I  GO  TO  20 

REAOI  07  .*  |   TPHOdF  (  I  I 

READ!  07  ,S  I   TNOZZ  I  I  I 

N  I  TEUS  =  1 
1  =  1*1 

go    ro    to 

20     CLiNTINUt 

*  w  I IE( h.iUI 

DO     JO      I = I .N I TEMS 

««[TeiO.50l  I.TAMBI  M.rPOOHEl  II  .TNDZZIII 
30     CONTINUE 

40     forma  r  I  tJX  .  •  AM8I  EN  T      I  C  I  •  ,  2  X  .  •  P  H  O  BE      (  C  I  '.  2X  .  'NOZZLE     ICI'I 
SO      FUMMA  I  (    1  X  .   I  J  .  OX  .2  |F6  .  J.oX   |   I 

srnp 

ENO 
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2.       TFIT 


tfi  t 


purpose : 


10 


21 


THIS     PROGRAM     PERFORMS      A     FIRST     URDEW     CURVE     FIT     BY      THE 
LEAST     SQUARES     METHOD     FOR     r  «  C     CCLumnS     OF     DATA.         THE 
FIRST     COLUMN     LISTS     VALUES     OF     X      anO     TmE     SECOND    LISTS 
VALUES     OF      F(X|      OR     Y.         FOR      THERMOCOUPLE     C»LI«R»I!UN. 
LET     X     =      THt      AMBItNT      TF "RER A r LRt .      ANO     LET     v     =     E I ThEH 

THt     PROBE     OR     NQ^Lt      TEMPERATLRE     aS     DESIRED. 


D  IMENSI 
SUMx=0. 
SUMxSU= 
Sum  xy  =  o 

SUMV=0. 

SU"NUM= 

SUMQEN= 

READ!  7  . 

DO     I  0      1 

REAO 

XSQI 

XV(    ( 

SUMX 

SUMX 

SUMX 

SUMY 

CON  T  I  NIJ 

xi  rt«b  = 

A=  (  SUMY 

•j-ixi  rc 

YtJAwr  bU 

DU     20      I 

YES  T 

*wl   IF  (  ft 

SUMM 

SUMO 

CG'U  I  NU 

RSQsSUM 

«K  I  TE  I  ft 

s  top 

End 


ON 
0 

0.  0 
.0 
0 
0  .0 

0  .0 
S  IN 

=  I  . 

(  7. 

1  I  = 

I  =  x 
=  su 

so  = 

Y  =  S 

=  su 

t 

FLO 
SSU 
MS* 

M  V  / 
-  I  . 
I  I  I 
.«l 

U«  = 

E  N  = 

t 

NUM 
ft  I 


x  I  150  1  .VI  1501  .XSOI  I  SO  I  . XY (  I  50  I  .  VEST (  ISOI 


I  TEMS 

N  I  TEMS 

*  )  X  I  [  I  .  V(  I  I 

x(  I  1**2 
I  I  ) *v (  I  ) 
MX ♦ X  I  1  I 
SUMXSU*  XSG  (  I  I 
UMX  V  ►  X Y  (  I  | 
MY*  V  I  I  | 

ATI  NI  TEMS  » 

MXSC-SUMX*SLMXY |/( X |TEMS*SUMXS0-SUMX**2( 

SIJMXY-SUMXSSUMV  |/(  X  ITEMS*SUMXOO-SUMX*^->  ' 

x  i  rt"S 

N  I  T£"S 

=  d*x I   1  )  ♦  A 

VtS T  (    I   I   .      V  (    I    I 

SUM  now*  (VEST  I  1  l-YH»w|«2 

SUMUEN*(Y( I l-VRAR|**2 


:2  I 


/SUMOEM 
A  ,i3  .3  SO 
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3 .       JETCURV 


c 

C  JETCURV 

C 

c 

C  PURPOSE:         THIS    PPOGRAM     FITS     JET     TRAJECTORY     DAT*     TQ     THE     EOUATION 

C  Z=ASY/(U»Yt      8V      THE     LEAST      SOUAWES     METHOO      AND     COMMUTES     A 

C  CORwOLATION     COEFFICIENT      IPJS2I.         DATA     POINTS      ARE     ftEAQ 

C 

C 

C 


THIS    PROGRAM     FITS     JET     TRAjECTCRV     DATA     TQ     THE     EOUATIOf 
Z=ASY/(U»Y|      8V      THE     LEAST      SOUAHES     METHOO      ANO     COMPUTES 
CORwOLATION     COEFFICIENT      IBJJ2I.         DATA     POINTS      ARE      REAC 
IN     FREE     FORMAT     FRQM     UNIT      "7"     CISK. 


O I  MENS  ION     Y(50l  .ZC50I.YV(50I.ZZ(50I  .YSOtSOl.YZISOl  .  ZEST (50  ) 

SUMY=O.U 

SUMYSO=0.0 

SUMYZ  =  0.0 

SUMZ=0.0 

SUMNUM=0 .0 

SUMUEN=0 .0 

REALX  7,t|NI  TEMS.FACTCR 

DO      10      1=1 .NITEMS 

REAOI  7.« I Y (  I   I . Z (  I  I 

Y (   I  I  =Y(   I   |/F  4CT0R 

YV( I |s| . 0/Y ( I ) 

ZZI  I  I = I .0/ ( Z ( I  l/FACTORI 

YSOI  I  ) =yy (  I  »s*2 

YZ I   I   ) =YY (   I  | *ZZ  (  I   I 

SUMY=SUMY*YY  (   I    | 

sumyso  =  sumysq ♦  ysu I  I  I 

SUMYZ  =SUMYZ»YZ  I  I  I 

SUMZ  =  SUMZ  *ZZ  (  I  I 
10  CONTINUE 

XlTc'«S  =  FLOAT(NI  TEMS  I 

AA=  I  SU^ZS'jUMYSQ-SUMYSSUMYZ  )/  (  X  I  TEMSSSL«YS0-SUMYC:=:2  I 
Ob=  (  X  I  TtMS*SUMYZ-SUMYC:SuMZ|/(  X  t  TE"S*S  t«VSQ-SUMYC:=2  I 
A=  |    . 0/  A  A 

e= A*fla 

ZHAR=^UMZ/X I TE"S 
OO     2  0      1=1. NITEMS 

ZES  r   (    I   I  stfOSVY  (    I    )  *  AA 

wR|Tt(o.J0IZZ(   II.  ZEST!   (I 

SUMNUM  =  SUMNUM»(ZfcSr  (  I  |-7BAR|J:S2 

SUMnEN  =  SUMOEN»  (ZZ1  I  l-«4HIS*2 

zo  con  r i Nut 

R  S  O  =  S  U  M  N  U  M  /  S  U  M  O  E  N 

1« I rE(6,40|4,U.hSC 
JO     FUTMAfllX.'Z      (ACTUAL!      =  •  .F  I  2. 8 . JX  .  •  Z      (EST)      :'.F12.fll 
*0     FORMA  1  (  1  X  .•  A=      •  ,F12.H,2X.'1:      •  .  F  1  2  .  a  .  2  X  .  «  R  **  ?.-      •   .F  1  2  .  6  I 

S  [HP 
END 
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4.       GRAB 


c 

C  GRAB 

C 

c  purpose:      oata    transfer    from    the    MP-<?e2eM  i  crocomputor 

C        TO  THE   IBM. 

C 

C        BY  BILL  CULtiRfcTH 

C        FOR  *E2*10.  FALL  QUARTER.   1982 

C 

C  FILEOfcF     05     TERMINAL 

C        FILtOEF   06   TERMINAL 

C        FILEOfcF  07  DISK  MVOATA  0  A T A   ( PERM  I 

C 

C        GLOBAL  TXTLIB  FORTM002  MOC2EEh 

C 

C  TYPE   tN  THE  AfiCVE  <*     LINES  TO  MAKE   THIS 

C  FORTRAN  PWOGRAM  RUN. 

C 

C 

DIMENSION  OATA(J00OI 

1  =  I 

"R I  TF (  6. 80  I 
80       FQRMAI  |  ^X  .  'QEG  IN   1NPUMNG  OATA  FRQM   ThE  MP-9826M 
C 

c 

10      CONTINUE 

REAUl  5. *l   OAT  A  I  I  I 

I  -  I  ♦  I 

IFIUMAI  I-U.NE.-20QI  GOTO  10 
C 

N| TfcMS  -   I- I 
h  FQRM A r  ( 2X ,  IS.  •    OATA  POINTS   *ERE  ENTERED. •! 

«H| IEI6.6I   N I TEMS 
C 

C  NOV*      THAT      ALL     DATA     HAS     4EEN     ENTEREO.      WRITE      IT     OUT     ON 

C  OlS*  . 

C 
<S  FO^MA  T  I  2X  .  »OA  T  A  (   •  .  I  5  .  •   I      =      '.IFIS.SI 

1  =  1 
20  «fi I  TE  I  7 .*  I      OAT  A(   f   | 

I  -  I  ♦! 

IF ( DA T A( I-| | .NF.-20C )      GOTC     20 
C 

C  ALL     UATA     HAS     HEEN      nRlTTFN     ONTO     DISK. 

C 

STOP 

ENO 
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5 .       TDATA 


c 

c 

c  ro* r » 

c 

c  purpose:       this   p«cc 

C  FCB»AT. 

c  tq   «e  r«  ic 

c 


Saw     PLACES      ri~e     aOCYANT     JET     3ATA     SECEIVEO 

-=-;32e    «  t  ccccp^  f?  -    ihto    a    »uoe    opofsl* 

IT  Ac5C  CONVERTS  UMTS  QF  LEsGT-  FROM  St 

0["ENSICN  XI2O0I  ,YI20OI  .ZI200I. TPOObE  (200)  .TAM6I200I  .TSOZ2I230I 
ISTDEVI200I  .*M»liOOI  . YM*(  200  I  ,Z«"I  2001 
1  =  I 
10  ^EiDI J7.SI   All) 

!FU(|  I.EQ.-200.0)   GC  TO  20 

««u  (  I  I =25  .4SX  (  r  I 

3t»U( U7,i|   Y ( [ ) 

Y««M  I  )  =25 .«  =  * I  I  ) 

READ! 07.il  lit) 

J"--  (  [  1=^5  .»=  *  (  t  I 

-- » ^i o? . *  i  rpnoee i t  i 

0£AJ(  G7,i|   TA-H  I  I  ) 

se»oio7,;i  r^cz^l  i  > 

^EiuU7.5|   STCEv  I  [  » 
s  i  TE«S=  I 
1  =  t  ♦  I 
GO  TQ   10 
20  CON  T  [  SUE 

»M|  TE (  6 . 40 ) 

Q ~   JO   1=1 . N 1 TE"S 

■airetb.sai  I  .x-«(  i  i  .  y«m(  t  ).Z"«(  r  i.tpw^e:  i  i,u"»i  t  i  .tnozzI  n 

1STDEVI   t I 

JO      CQNtlNUt 

*0     FQt»«  a  r  (  a*  .  •  x     |»»l',4«,  'y     [am|i.AX.«2      (  w»  |  •  ,  2  x  ,  ■  poo3E      ICI'.l*. 

I'AMblENI      ICI"  ill.'NCZOc      ICI" .IX. «ST0      ~cv      !C>*' 
50     cG4v.ar(1x.tJ.JX.2|F7.j.jx).>--7.3.4x,..(F*,.j.*xi| 

srnp 
£■•: 
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C0NT0UR4 


c 

C  CONTOUR* 

C 

C 

C  PUHPOSe:         THIS    PROGRAM     IS    OFSIGNED     re     CISPL*Y     8uor»Nr     jFT 

C  TEMPfcRATuRE     DATA     USING     A     CONTOUR     PLOTTING     PACKAGE 

C  AVAIL^eLE      alTH     OIS5PL*.         Tt-F     FIjLLOwInG     RAW     OATA      IS 

C  RFAO     FRCM     A     O  I  S*     PILE:     PWOEF     POSITIONS     IN     XYZ 

C  COORDINATES     RELATIVE      TQ      THE      NOZZLE      TIP.      AMBIENT 

C  TEMPERATURE.      PROBE      TEMPERATURE      ANO     NQZZLt      IFmPFRATUPE. 

C  THt"     PROPE      ANO     NOZZLE     TEMPERA  ruwES     ARfc      TRANSFORMED 

C  BY     CALIBRATION     COEFFICIENTS      ANO     NORMALIZEO     WW T      AMBIENT 

C  TEMPERATURE.         THE      XYZ     COURCINATES     ARE     CONVERTED      TO 

C  XSw     COORDINATES     RELATIVE      TC      IHt     INTERSECTION     OF      THE 

C  UA  r A     PLANE     «ITh     THE     CENTERLINE      TRAJECTORY     UF      THF      JET 

C  (THE     S-AxIS      IS     TANGFNT      TO      THt      TRAJECTORY!.          THIS 

C  PRUGRAM     ALSO     COMPUTES      THt      SAIF     OF      mFaT      TRANSFER     FROM 

C  THE     j£T      TC     THE      AMBIENT. 

C 

C  NOTE:          THE     FOLLOWING     VALUES     MUST     BE      INSERTED     AS      THE 

C  FIRST     LINE     OF     CATA      IN     FREE     FCR-AT:      THE      TOTAL     NUMBER 

C  OF     DATA     POINTS     (NITEMSI.      Tt-F     ACUTE     ANGLE     BETWEEN 

C  THE     OATA      PLANE      AND     HORIZONTAL      I ThE T A |      ANO      THE      VERTICAL 

C  OISTANCE     BETwdfcN      Thf      Z-AXIS      ANO      THE      INTERSECTION 

C  POINT     CISCUSSEC     ABOVE      (ZAI,      THE     CENTERLINE     VELOCITY 

C  IN     M/s      (VELl      *NO      JET     WIDTH      IN     MM      (WIDTH). 

C 

c 

OIMENSION     X  (  100  I  .  Y I   1001.  TIIOOI.TP  (1001.  TA(   100  I  .  TNI   100I.TMATI   10.101 

COMMON   wORKIlbOOOl 
C 
C 

C*  *  *  O  READ  OATA  FROM  UNIT   ■•  T "  O  I  SK  *   *  S  * 
C 
C 

REAOI 7.S  I    NlTEMS.THETA.ZA.vEL.wlCTH 


00     20      1  =  1  .N  I  TEMS 

HEA0|7.*|  i,X|  I  |.A,Y[  II.TPIII.TAI    I  I   .TN(  I   I   .  A 

Y(l|       =       (  Y  (   I    |-ZA   |  /SI  M    THET A( 
C 
C  CALIBRATION     COEFFICIENTS 

TP(   I  )  =  .9B4e0Qe75*TP(  I  )  ♦  I  .7*U7<)bC5 

TN(  I  l=.9  13223*2b*TN(  [  |*|.d*75»ne6 
C 

Till   =   (TP(|  I   -   TA  I  I  I  I  /(  TNI  I  I   -   TA||  II 
20    CONTINUE 
C 

C   30    FORMAT  (  2X  .  •  I  .  X.  Y.  T :  •  .  15 .  3F  I  5.5  I 
C 
C 

C*  *  *  *  DIMENSION  OF   TMAT   S  S  C  * 
C 

c 

I XO I M   =   10 
I  YD  I M   :   10 

c 
c 

c 

CS   *   S   4!   FIND   THF   MAXIMUM   ANO   MINIMUM   TEMPFftATURES   -   S   *   S 

c 

c 

TMIN     =      I  00  .0 
TMAX       =       0.0 
XM A  X=  0  .  0 
YM A  X  =  0  .  0 
XM(N=I00.0 
YM  I  N=   I  00  .  0 

c 

00     40      I = 1 .N I TfcMS 

[F|T|I|        .GT.       TMAXl  TMAX      =      T(I) 

iF(r(u     .Lf.     r  x  i  n  i        r»iN    -     tiii 
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c 
c* 

c 

c 
c 
c* 

c 
c 


IF  (  X  (  I  I  .GT.Xy4X|X««X:X||  | 
IFIY|  I  I  ,(,r.VM«X|yW»X:»|  I  ) 
IF|X||I.LT.XN|MXMIN:X|1| 

IF  (  Y  (   I    I    .LT  ,rMlM»M|N:r|   [   I 
40        CQNMNUc 

:  *  UETERWINt"   THE   MAXIMUM  VALUES  OF   X  AND  Y  «(TH1N  A 

OUAQRANT.   ASSUMING  ThE  STREamwISE  AXIS   IS  CENTERED 
IN  A  SUUARE  MATRIX  PLANE 

0X=(  [XMAX-XMINl/FLOAri  IXO  l«4  )  1/1000. 
OY=(  (YMAX-YMIN|/FLQAT(  I  YO  [M  |  I  /  I  00  0  . 

OuAOX=|FLUATI IXUIM|/2.Q|*CX 
OUAUY=|FLOAT( IY0IM|/2.0)*CY 

!   «  F  I  NO  THE   INCREMENTAL   AREA  FOR  HEAT   TRANSFER  CALCULATIONS 

AREA  =OX*0 Y 

*  *     S     ESTABLISH     rue.      INCREMENT     SIZE    FQR     CONTOUR     PLOTS     3     S     S     * 
TINCR=( TMAX-TMINI/5.25 

*  *     «     GENERATE     THE     GRID     *    *     *    * 


C 

C* 

C 

C 

c 


c 

c*    * 

c 


je  r*» 

2  ..<*  I 


CALL  COMORS 

CALL  TEk.6  1  e 

Call  PAGE  I     d.  . fl .  I 

CALL  BLU«UP      (0.75  1 

CALL  PHYSOR ( I .25. 1 . I 

CALL  APEA2015.5  .5.5  I 

CALL     HE 1GHT  ( 0.  1 00  I 

CALL     C»RTUG 

CALL      INTAXS 

CALL      XNAMfc(.X      (MM|*«,100> 

CALL      YNAMEI'W       (MM)t«,l00| 
C  ALL     XT  ICKSI  2  I 
C  ALL     H  I  C*  SI  2  I 

CALL       Y  AX  ANG  (   <50  .  0  ) 


CALL  SwISSM 

CALL  HE AO I N ( • TEMPER A TLRE     CONTOURS      IN      A     fcUUYANT 

CrtLL  mFA01N(»«ITH      A     CROSSFL O* INC     AMH  I  EN  1% •   .  1  00 , 

CALL  hEaOIN('(25X     FlC«     R  «  IE  I  *  •  .  I  0  C  .  I  .  .  4  I 

CALL  HEAD  IN (  •  (PL ANE  A  |  s •  .  I  00  .  I  .  .4  | 

CALL  GRAF ( -JO  . 0.05.0. JO. 0. -JO. 0.05. 0.20.  01 

it     U     $■     GENERATE      THE      INTERPOLATED      TEMPfcRATLPt     MATRIX     "fAT1 

CALL     dCOMUN ( I oOOO I 
CALL     Ze«Sfc ( TMIN I 

CALL     PGNMA  T  (   1 XO I M  ,  I  YD  I  M  ) 

CALL      GETMATIX.  Y.T.N  ITEMS. 01 

CALL     F NOMA  T I   TM A T  , 0  I 

*     *     *     SOLVE      THt     RATE     OF     HEAT      TRANSFER     S     *      *     « 

USUM=0 . 0 

UO      70      I=l.lXUIM 

00     SO     J: I   .  I  YO  I  M 

RO=  (   I  .  0000*  J  J     »     3.5  I  9  196-SS  r»A  T  (   [ 
1  -*>.9;jft*:oa7E-eSTMAT|t.J|SS2 

1  »      J.4l2')o?qE-PiTM4|(   l,j|-*J 

I  -'J.7<»Os7tF-tlST-4AT(l.J)SS<»|StOOO. 

XX-QUAUX-FLOATI  J-l  I  "SOX 
Y  Y=i)o  ai;y-FlO  a  T  (  I-i  |~ly 
R AO luS -SLR T ( x XCi 2      ♦      YYSS2I 


100  . 2 .  . 4  I 


J  I 
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60 
/0 


ao 


V=VELSEXP(-Ra0  1US*S2/v«IDTHS«2) 

cont i Nut 
cont i nue 
*R i re ( G.ao iosum 

FOUmaT l IX. «0=   '.F20.5I 


c 

c* 

c 


PLOT   The  CONTOURS  S  *  *  * 


CALL 
CALL 
CALL 
CALL 

Call 

CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
STOP 
END 


CONMAK|TM«r. I  XO  t  *, 

CONL IN| 0 , • SOL ! US' , 

CON anG  (  90  .  O  I 

KASPLNI  0 . 2b  I 

CONTUM ( 1 . 'LABELS • , 

ME  ICriT   (  0  .  05  I 

WLMESS ( 'EXP   | C  t • , 

DOT 

CMIOI I .11 

tNOPL  (  0  t 

OONEPL 


I  YOIM.  Tf NCR  » 
•LARELS* .1 .101 


•ORA\»  •  ) 

I  00 .-59, -59  I 
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APPENDIX  D 
TABULATED  DATA 
TABLE  1 
ROTOMETER  CALIBRATION 

%  Flow  ml/s  Std  Dev  (ml/s) 

10 
15 
20 
25 
30 
35 
40 
45 
50 
55 
60 
65 
70 
75 


5.94 

0.32 

8.02 

0.56 

10.05 

0.41 

11.85 

0.24 

14.03 

0.33 

15.98 

0.58 

17.41 

0.31 

19.09 

0.34 

20.75 

0.23 

22.38 

0.28 

24.27 

0.22 

26.03 

0.20 

27.64 

0.37 

29.53 

0.24 
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TABLE  2 
TEST  RESULTS 


Crossflow  velocity: 

Nozzle  flow  rate: 

Nozzle  inside  diameter: 

Nozzle  discharge  velocity  (mean) 

Nozzle  temperature  (mean) : 

Ambient  temperature  (mean) : 

Froude  number: 

Michaelis-Menter  equation: 
Coefficient  A: 
Coefficient  B: 


.13  m/s 
11.85  ml/s  (25%) 

7 . 144  mm 
29.558  mm/s 

41.8  °C 

24.9  °C 
14.8 

2.64284325 
1.03698254 


Plane 

Y  (mm) 

* 

(degrees) 

Q  (W) 

A 

7.327 

46 

1.20 

B 

21.370 

58 

6.392 

C 

39.688 

70 

15.349 

D 

62.889 

80 

18.388 

E 

87.313 

86 

27.628 
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